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Abstract 


The perturbation of the lunar motion caused by a hypothetical violation of 
the equivalence principle is analytically worked out in terms of power series a 
la Hill-Brown. The interaction with the quadrupolar tide is found to amplify 
the leading order term in the synodic range oscillation by more than 62 %. 
Confirming a recent suggestion of Nordtvedt, we find that this amplification 
has a pole singularity for an orbit beyond the lunar orbit. This singularity is 
shown to correspond to the critical prograde orbit beyond which, as found by 
Henon, Hill’s periodic orbit becomes exponentially unstable. It is suggested 
that ranging between prograde and retrograde orbits around outer planets 
might provide future high precision orbital tests of the equivalence principle. 
It is argued that, within the context of string-derived non-Einsteinian theories, 
the theoretical significance of orbital tests of the universality of free fall is to 
measure the basic coupling strength of some scalar field through composition- 
dependent effects. Present Lunar Laser Ranging data yield the value 7 = 
(—1.2± 1.7) X 10“'^ for the effective Eddington parameter 7 = 7 — ! measuring 
this coupling strength. 
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I. INTRODUCTION 


Gravity seems to enjoy a remarkable universality property: all bodies are experimentally 
found to fall with the same acceleration in an external gravitational held, independently of 
their mass and composition. Although Galileo |l| was the hrst 0] to suggest in a clear and 
general way that this property of universality of free fall might hold true, it was left to Newton 
0 to realize the remarkable conceptual status of this universality: exact proportionality 
between a particular force (the weight) and the general dynamical measure of inertia (the 
mass). Newton went further in performing the hrst precise laboratory tests of the universality 
of free fall (pendulum experiments; precision ~ 10“^). It is less well known that Newton 
went even further and suggested to test the universality of free fall of celestial bodies by 
looking for a possible miscentering of the orbits of satellites around Jupiter, Saturn and the 
Earth |^. More precisely, Newton considers a possible violation of the ratio weight {w) over 
mass (m), i.e. 


^12 = 


{w/m)] 


{w/m)r 


1 ^ 0 , 


( 1 . 1 ) 


where 1 labels a satellite and 2 a planet, the weights wi, W 2 being the gravitational forces 
exerted by the Sun (label 3). He says, without giving any details, that he has found “by 
some computations” that the center of the orbit of the satellite 1 around the planet 2 will 
be displaced (in the Sun-planet direction and away from the Sun if S 12 > 0) by the amount 


(5|Xc - X 3 I = +-Svia' , 


( 1 . 2 ) 


where a' denotes the radius of the orbit of the planet 2 around the Sun. In modern phrase¬ 
ology, one can say that Newton predicted a “polarization” of the satellite’s orbit in the 
Sun-planet direction (away from the Sun if 812 > 0). Then Newton used his theoretical 
estimate (|1.2|) to conclude from the observed good centering of the orbits of the satellites of 
Jupiter that |(5i2| < 10“^, a number comparable to the result of his pendulum experiments. 
Actually, this upper limit obtained from Jovian satellites is wrong, as Newton’s theoretical 
estimate (|1.2|) is incorrect both in magnitude (being a gross overestimate in general) and in 
sign (see below). We could not End any information about Newton’s original calculations in 
his published papers. It is surprising that Newton did not remark that, as a consequence of 
his estimate (O), a value IJ 12 I = 10“^ would also entail an unacceptably large polarization 
(one fifth) of the Moon’s orbit. 

As far as we are aware, Laplace was the first to realize that the best celestial system 
to test a possible violation of the universality of free fall ( |1.1|) is the Earth-Moon system 
(1 = Moon, 2 = Earth). In P| he derived a rough estimate of the main observable effect of 
S 12 on the angular motion of the Moon. Then, he noticed that even a very small ^12 ^ 0 
would spoil the agreement linking his theoretical derivation of the “parallactic inequality” of 
the Moon, the set of observations of the lunar motion, and the direct measurements of the 
solar parallax. He concluded that an upper bound to the fractional difference in acceleration 
of the Moon and the Earth toward the Sun is 


|<^ 12 | < 


3410000 


~ 2.9 X 10 


-7 


(1,3) 
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It is remarkable that the limit ( p..3|) is much better than the limit (I^abI < 2 x 10" 
obtained some years later by Bessel through improved pendulum experiments @], and has 
been superseded (though not by much) only by the work of Eotvos in 1890 (|^ab| < 5 x 10“®; 
0). [The later results of Eotvos, Pekar and Fekete improved the bound to 3 x 10“®; [|].] 
As we discuss in Appendix D, in spite of some obscurities in his reasonings and the lack of a 


fully accurate calculation of the effect of S 12 in longitude, Laplace’s hnal bound (|1.3|) turns 
out to be a conservative upper limit, given the information he had. 

In 1907, Einstein 0 deepened the conceptual implications of the property of universality 
of free fall by raising it to the level of a “hypothesis of complete physical equivalence” between 
a gravitational held and an accelerated system of reference. This heuristic hypothesis was 
used very successfully by Einstein in his construction of the theory of general relativity, and 
became later enshrined in the name “principle of equivalence”. 

Within the context of relativistic gravitational theories, the use of the Moon as a sen¬ 
sitive probe of a possible violation of the equivalence principle for massive, self-gravitating 
bodies has been rediscovered by Nordtvedt in 1968 [T^. His idea was that self-gravitational 
energies might couple non universally to an external gravitational held in theories having a 
diherent structure than general relativity |n[, |^. [Let us note that Dicke had mentioned 


this possibility earlier [^, |^, [^ but had not explored its consequences in detail.] Any¬ 
way, Nordtvedt, unaware both of the old ideas of Newton and Laplace, and of the more 
recent ones of Dicke, realized that the planned Lunar Laser Ranging (LLR) experiment was 
providing an exquisitely sensitive tool for testing the universality of free fall of massive bod- 
T0[| . Performing a /irsf-order perturbation analysis of the lunar orbit (assumed circular 


les 


and planar) in presence of a violation of the equivalence principle, 612 7 ^ 0 (see Eq. 0 
with the labels 1 and 2 denoting the Moon and the Earth, respectively), he provided the 
hrst analytical estimate of the corresponding range oscillation: 


(dr)^^^ = C‘^^^ 5 i 2 a' cos [{n — n')t + tq] 


(1.4) 


with 


^( 1 ) _ 1 + 2n/(n-n0 2 

n2-(n-n')2 


(1.5) 


Here, n denotes the (mean) sidereal angular velocity of the Moon around the Earth, n! the 
(mean) sidereal angular velocity of the Earth around the Sun, and a! denotes the radius of 
the orbit of the Earth around the Sun (assumed circular). The angle r = {n—n')t+TQ is equal 
to the difference between the mean longitude of the Moon and the mean longitude of the 
Sun (as seen from the Earth). For completeness we derive Eqs. ( p..4| ), (|1.5| ) in Appendix A. 

Let us note in passing that the cos r dependence of the range oscillation ( [1.4|) is equivalent 
(when disregarding the perturbations of the motion in longitude) to displacing the center 
Xc of a circular lunar orbit in the Earth-Sun direction (toward the Sun if 612 > 0) by the 
amount 


(5|xc - X 3 I = 


12® 


( 1 . 6 ) 


The result ( [1.2|) of Newton can therefore be viewed as implying a range oscillation of 
the type (|1.4|) with (^Newton = — 1 independently of n and n' . By contrast, the first order 
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estimate (|1.5| ) contains the small dimensionless parameter[| 


n 


m = 


n — n' 


(1.7) 


which is m ~ 1/12.3687 for the Moon and much smaller for the (Galilean) satellites of 
Jupiter (e.g. m ~ 3.86 x 10“^ for Jupiter IV). More precisely, Eq. (ED can be rewritten as 


(1) o 1 + 3m 

’ = -m -T— = -m 


3 

—'i 

2 


1 + fm 


3 

—t 

2 


. 1 1 2 

1 H —m - m 

6 12 


+ 


( 1 . 8 ) 


In 1973, Nordtvedt suggested that a more accurate value of the coefficient C in the 
COST (or “synodic”) range oscillation (1^) would be obtained by replacing, in the denom¬ 
inator of Eq. (|1.5|) , the angular velocity n by the frequency of radial perturbations: 

0 


^rad = I = cn where the perturbation series giving c reads (see e.g. 


0 ) 


c = 1 




= 1 


-m 


n 


177 


-m 


(1.9) 


[In the case of the Moon, the series (|1.9|) is very slowly convergent. The full value of 
1-c ~ 0.008572573 13 is more than twice the lowest-order correction |m^.] The correction 
of Ref. amounts numerically to increasing the hrst-order result (|1.5|) by about 13 %. 

In 1981, Will [|^] tried, more systematically, to estimate the higher-order corrections in 
the coefficient C due to the mixing between the perturbation (1^) at frequency n — n' and 
the tidal perturbations at frequencies 0 and 2(n — n'). He suggested that the hrst order 
result should be multiplied by a factor 1 -|- 2n'/n = 1 -|- 2m -|- Oijn?), i.e. amplihed by about 
15 %. As a result of these (coincidentally equivalent) prescriptions, the literature on the 
0, 0, 


“Nordtvedt effect ” 0, 0, 0 has, for many years, used as standard estimate 

for the range oscillation 6r = cost a value C ~ 1.14 [corresponding to about 

9.377COST meters in metrically-coupled theories; see below]. 

Actually, both modihcations (suggested in 


and 


of the hrst-order result fall short 


of giving an accurate estimate of the effects due to higher powers of m. In fact, they do not 
even give correctly the second order in m. For completeness, we compute in Appendix A, by 
the standard perturbation theory of de Pontecoulant the contribution at order 0{m?) 


and hnd that it amounts to multiplying the hrst-order result by 1 -|- |m + 0{m 


i.e. 


(^( 1 ) ^ (j{2) ^ 


- 


1 + 

-m + 0 (m^) 


2 ^ 

3 

r. 14 

—m 

1 -b —m -b 0 1 

2 

L 3 


C'(i) 


( 1 . 10 ) 


By contrast, [0 and []0| would give 


and 


12 


6 ’ 


respectively, for the coefficient of m in 

14 „ 


the correcting factor within square brackets in the second Eq. (|1.10|) . Note that ym = 37.7 


^ This is the version of the small parameter which is appropriate to our Hill-Brown treatment. 
Beware of the fact that the more traditional perturbation approaches denote by the letter m the 
quantity m = n'/n. 
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% for the Moon. Recently, Nordtvedt argued that the full perturbation series in powers 
of m (created by the interaction with the orbit’s tidal deformation) will cause a rather large 
numerical amplihcation of the synodic oscillation (|1.4|) . However, he gave only an incomplete 
theoretical analysis of this amplihcation. The only literal result he gives (his Eq. (2.33)) is 
equivalent to the second-order effect derived in Appendix A and leading to Eq. (|1.10|) . We 
could not extract from his paper a precise value for the full effect of the m-series. 

The aim of the present paper is to provide, for the hrst time, a full-hedged Hill-Brown 
analytical treatment of the orbital perturbations caused by a violation of the equivalence 
principle. Our results will allow us notably to give a precise numerical value for the full 
range oscillation in the case of the actual Moon^. Namely, we obtain below 


5r = 2.9427 x 10 ^^ 5 i 2 cosr 


cm 


( 1 . 11 ) 


corresponding to a full coefficient C = |m x 1.62201 which is larger than the hrst-order 
value (|1.8|) by more than 60%. More generally, we shall be able to discuss in detail the 
dependence on m of the range oscillation: see Eqs. (p.60|) - (|2.62|) and Appendix B. These 
results are summarized in Fig.l below. Our results conhrm the suggestion of Ref. that 


when m increases (corresponding to prograde orbits beyond the actual lunar orbit) the cost 
range oscillation eventually becomes resonant and is (formally) inhnitely amplihed. We have 
some doubts, however, about the practical utility of such a resonant orbit, notably because 
we show that it occurs precisely at the value 


m = rucr-= 0.1951039966... (1.12) 

(corresponding to a sidereal period Ter = '^cr/(l + ^cr) yr = 1.95903 month) where the 
orbit becomes exponentially unstable. Armed with our theoretical understanding of the 
m-dependence of the cost oscillation, we suggest other orbits that might be practically 
interesting (retrograde orbits, orbits around outer planets). Finally, we emphasize that 
within the context of modern unihed theories the most probable theoretical signihcance of 
orbital tests of the universality of free fall is the same as that of laboratory tests, namely to 
measure, through composition-dependent effects, the strength of the coupling to matter of 
some long range scalar £eld(s). The basic measure of this coupling strength is embodied in 
an effective Eddington parameter 7 = 7 — ! which governs both the standard post-Newtonian 
effects (including the violation of the strong equivalence principle cx: 7 = 4/? — 7 ) and the 
composition-dependent couplings (violation of the weak equivalence principle). Actually, 
string theory suggests that the former contribution (proportional to the gravitational binding 
energy) is negligible compared to the one due to a violation of the weak equivalence principle. 
Interpreting the latest LLR observational results within a recently studied class of string- 
derived theoretical models, we conclude that present orbital tests give the excellent constraint 
7 = (—1.2 ± 1.7) X 10“^. [This limit is comparable to the (similarly interpreted) constraint 
coming from laboratory tests of the weak equivalence principle: 7 = (—0.8 ± 1.0) x 10“^.] 


^Note, however, that we consider only the Main Lunar Problem, i.e. that we neglect the terms 
proportional to the squares of the lunar and solar eccentricities, and to the square of the lunar 
inclination, which are expected to modify our numerical estimates by % 1 %. 
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The plan of this paper is as follows. Section II presents onr Hill-Brown approach and gives 
the analytical resnlts obtained with it. Section III discnsses the physical conseqnences of onr 
resnlts. Many technical details are relegated to some Appendices: Appendix A presents the 
standard de Ponteconlant treatment of Innar theory and nses it to derive the second-order 
resnlt ( |1.1CI| ), Appendix B gives some details of onr Hill-Brown treatment, Appendix C treats 
the link between certain commensnrabilities of frequencies, linear instability and the pres¬ 
ence of pole singularities in perturbed motions, and hnally Appendix D discusses Laplace’s 
derivation of the remarkably good limit (|1.3|) on Eu- 


II. HILL-BROWN TREATMENT OF EQUIVALENCE-PRINCIPLE-VIOLATION 

EFFECTS 

A. Introduction 


Relativistic effects in the lunar motion have been investigated by many authors. The 
pioneers in this held are de Sitter (who computed the general relativistic contributions to 
the secular motions of the lunar perigee and node as observed in a global, barycentric frame) 
and Brumberg ||^ (who gave a comprehensive Hill-Brown treatment of the post-Newtonian 
3-body problem). Later works studied non-Einsteinian effects, notably those associated with 
the Eddington post-Newtonian parameters (3 and 7 . The most comprehensive and accurate 
analytical study of post-Newtonian effects in the lunar motion (described in a barycentric 
frame) is due to Brumberg and Ivanova |^. For general accounts and more references see 
the books and [^. Let us also mention the semi-analytical treatment of the general 


relativistic perturbations of the Moon by Lestrade and Chapront-Touze pO 

0, m 


However, apart from the work of Nordtvedt 


the studies of non-Einsteinian 
effects in the lunar motion have not considered the effect of a violation of the equivalence 
principle. The results of the present paper can therefore be considered as a completion of 
Ref. which gave an accurate Hill-Brown theory of all the other Einsteinian and non- 
Einsteinian effects. In fact, as pointed out long ago by Nordtvedt, the effects of a violation 
of the universality of free fall are the most prominent non-Einsteinian effects in the lunar 
orbit, and therefore deserve an accurate study. Indeed, most of the non-Einsteinian effects 
are non-null effects, i.e. correspond to modihcations proportional to/3 = /3 — lor 7 = 7 — 1 
of observable relativistic effects (as seen in a local, geocentric frame) predicted by Einstein’s 
theory. As the latter are at the few centimeter level pT 


2 ^ 


, which is the precision 

of the LLR data, they can be of no use for measuring (3 and 7 at an interesting level (say 
< 10“^). An exception must be made for secular effects and for the parameters describing the 
temporal and spatial transformation linking a local, geocentric frame to a global, barycentric 
one, e.g. the parameters entering the de Sitter-Fokker (“geodetic”) precession. [Recent work 
concludes that geodetic precession alone constraints 7 at the 1 % level]. 

Besides the “Nordtvedt effect” proper (i.e. the effect of ^12 7 ^ 0), that we discuss here, 
there are some other null effects which are more sensitive to f3 and 7 than the non-Einsteinian 
modifications of non-null general relativistic effects. A subdominant null effect comes from 
the violation of the equivalence principle associated with the gravitational binding energy 
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of the Earth-Moon system. In lowest approximation (linear in m), it is equivalent (see e.g. 
Eq. (3.14b) of Ref. 0) to replacing ^12 by 

■ 5.2 = (“)' , ( 2 . 1 ) 

where a denotes the semi-major axis of the lunar orbit and where rj denotes, as usual, the 
combination 


77 = 4/3 — 7 = 4/3 — 7 — 3. 

In general, ^12 is the sum of two physically independent contributions 

( 7-igrav pigrav \ 

-S-^j ■ 

mic^ m2C^ J 


( 2 . 2 ) 


(2.3) 


The first contribution ^12 = 61 — 62 is generically expected to be present because the best 
motivated modihed theories of gravity violate the “weak equivalence principle”, i.e. contain, 
besides Einstein’s universal tensor interaction, some composition-dependent couplings that 

36|). The second contri- 


make laboratory bodies fall in a non universal way (see e.g. 
bution on the right-hand side of Eq. ( p.3|) (proportional to r] like the correction in Eq. (|2.1|) ) 
contains the gravitational self-energy of the bodies {A = 1,2) 


E^^'' = —{G/2) f f d^xd^x'p{^x)p(x')/\:>^ 
Ja Ja 


(2.4) 


and is the one hrst pointed out by Nordtvedt n. H- As indicated by Dicke |T^, it is 


present in all gravity theories where the effective, locally measured gravitational “constant” 
may vary from place to place (see e.g. Section V B of [^). We shall take as nominal 
values for the gravitational self-energies of the Moon (label 1) and the Earth (label 2) the 
values adopted by Williams, Newhall and Dickey [BSl, namely /mic^ = —0.19 x 10“^°, 
Er^lm 2 C^ = -4.64 x 10-^°, so that 


^grav pgrav 

-^ = 4.45 X 10 “^° . 

rriic^ m2C^ 


(2.5) 


We then End numerically that the modihcation due to the gravitational binding energy 
of the Earth-Moon system, —^pri^a^/c^ = —^pG{mi +m 2 )/ac^ in Eq. ( |2.1|) , is (to hrst 
order) equivalent to decreasing the nominal value (|2.5|) by —0.039 x 10“^°. This represents 
a fractional change of (|2.5|) by —0.87 % which is probably smaller than the uncertainty in 
the estimate (|2.5| ) associated with our imperfect knowledge of the internal structures of the 
Earth and the Moon. These orders of magnitude illustrate the fact that the overwhelmingly 
dominant sensitivity of the lunar motion to non-Einsteinian effects comes from the terms 
proportional to 612 that we concentrate upon in the following. 


B. Three-body Lagrangian 

The Lagrangian describing the Wbody problem in the currently best motivated rela¬ 
tivistic theories of gravity, i.e. those where gravity is mediated both by a tensor held and 
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a scalar field with, generically, composition-dependent couplings (see e.g. 
written as 


|]36|), can be 


-^/ 3 , 7,(5 — ^GR + Li3 + , 


13 -r 


( 2 . 6 ) 


where Lgr denotes the general relativistic contribution (in which one should use an ef¬ 
fective value of the gravitational coupling constant G which incorporates the composition- 
independent part of the interaction mediated by the scalar held), 


1 _ ^ GniAmB , .2 

^7 = ^7 E ’ 


(2.7) 


denotes the (non-tensor-like) velocity-dependent part of the two-body scalar interaction 
(one-scalaron exchange level). 


_ I - G^rriAmBmc 

^0 ~ 2 ^ 




tabTac 


( 2 . 8 ) 


denotes the modihcation of the non-linear three-body general relativistic interaction due to 
the scalar interaction, and 


L-s 


11: (s.,+4) 


GuIaITIb 

Tab 


(2.9) 


with 






^grav 


9 ’ 

rriAC 


( 2 . 10 ) 


represents (to lowest order) the combined effect of the composition-dependent couplings 
{5a 7 ^ 0; violation of the “weak equivalence principle”) and of the Dicke-Nordtvedt con¬ 
tribution due to the spatial variability of the effective gravitational coupling constant 
(?7 = 4/9 — 7 7 ^ 0; violation of the “strong equivalence principle”). For a direct, held-theory 
derivation of and and the expression of the phenomenological Eddington parameters 
7 = 7 — 1 and ^ = (5 — 1 in terms of the basic coupling parameters of the scalar held (as 
well as the generalization of these results to the case of strongly self-gravitating bodies) see 


Ref. 

We assume here that all the general relativistic contributions to the lunar motion (and 
to its observation through laser ranging) are separately worked out with sufficient accu¬ 
racy, using, for instance, the new, complete framework for relativistic celestial mechanics of 
Ref. (which provides the hrst consistent relativistic description of the multipole moments 
of extended bodies). Following the discussion above, we henceforth discard the subdomi¬ 
nant contributions coming from and to concentrate upon the ehects due to L^. [The 
barycentric frame contributions oi and have been accurately computed by Brumberg 
and Ivanova |T^ and can be linearly superposed with the ones of Lg.] Finally, it is enough 
to consider the sum of the lowest-order approximation to Lgr and of Lj, namely 









, , 1 2,1 Gab'i^a'i^b 

L{xa, ^a) = 2^ t^^a^a + o 


AjtB 


Tab 


( 2 . 11 ) 


where 


Gab = G 


1 + + ^. 


B 


( 2 . 12 ) 


with 5a of the form given in ( 2.1U| ), represents the effective, composition-dependent gravi¬ 
tational coupling between the massive bodies A and B. In Eq. ( p.ll| ) = d'^A/dt denotes 

the (barycentric) velocity of body Z, uia the (inertial) mass of A and tab = \^a — the 
distance between A and B. We consider a three-body problem, and more particularly the 
Moon-Earth-Sun system (1 = Moon, 2 = Earth, 3 = Sun). Evidently, all our analytical 
results will apply if 2 denotes another planet, and 1 one of its natural or artihcial satel¬ 
lites. [Note, however, that the relative orders of magnitude of the non-Einsteinian effects is 
different for low-orbit artihcial Earth satellites. See for a recent discussion.] 

Starting from Eq. (|2.11|) , we hrst separate the variables describing the motion of the 
center of mass of the Earth-Moon system. 


moXo = miXi -F m2X2 , 
mo = mi -t- m2 , 


( 2 . 13 ) 

( 2 . 14 ) 


and vq = dxo/dt, from those describing the relative lunar motion. 


Xi2 = Xi - X2 , 

mim2 

hl2 = -^- , 

mi -|- m2 


( 2 . 15 ) 

( 2 . 16 ) 


and Vi 2 = (ixi 2 /(it. This yields 


rf 2,1 2,1 2,^ ^irn2 mimg m2m3 

L(Xo, Xi2, Xs; Vo, Vi2, Vg) = -moVo -/U12V12 -mgVg -F Gi2 - h 613-h G23 - 

2 2 2 ri2 ri3 r23 

(2.17) 


where 


ri3 = |x3 - xo - X2X12I , ( 2 . 18 a) 

?^23 = |x 3 -xo-FXiXi 2 | , ( 2 . 18 b) 

Xi = mi/mo , X2 = m2/mo = 1 — Xi . ( 2 . 18 c) 

Expanding and r23 in powers of ri2/r3o (rgo = [xgol with xgo = xg — xq) leads to 

L = Lo3(xo, Xg, Vo, Vg) -|- /ii2Tl2(Xi2, V12, Xgo) , ( 2 . 19 ) 

LqS = T^^iQVo + Tt^^SVg -I- (XiGig -I- X2G23) -- - , ( 2 . 20 ) 

1 I ro 3 

Ti2 = 77V12 + Gi2- \- Ri + R2 + R3 + ■ ■ ■ , (2-21) 

2 ri2 
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( 2 . 22 a) 


Rl — m 3 (—Gi3 + 6 * 23 ) 3^12 - 

^^30 

R 2 = TTT^S (^13X2 + G23-^i) 3^12^12 - ) ( 2.22b) 

^"30 

Rl = bns (-GaXl + ChiXf) x\,xl,x';, 0®— , (2.22c) 

3! ^ ^ ^ ^30 


where = d/dx^, d^f = d’^/dx^dzi, _ Note that the suffices 1,2,3,... in R^s have 

nothing to do with the body labels A,B = 1, 2, 3, but keep track of the successive powers 
of Xi 2 . 

To a very good approximation we can consider that, in the (normalized) Earth-Moon 
Lagrangian (|2.21|) , the motion of the Sun with respect to the Earth-Moon barycenter, X 3 o(t), 
is obtained by solving the two-body Lagrangian (S)- After separating in Eq. ( p.20|) 
the motion of the center of mass of the Earth-Moon-Sun system, the reduced Lagrangian 
describing the dynamics of the relative motion X 30 is 


where we introduced 


f relative 
-^03 


1 2 , ^ mo + 1713 

2V„3 + G„3—— 


(2.23) 


Go 3 = X1G13 -|- X2G23 ■ 


(2.24) 


Therefore, seen from the Earth-Moon barycenter, the Sun undergoes a Keplerian, elliptic 
motion corresponding to a total effective gravitational mass Go 3 (^^o + ^ 3 ) = Go 3 (’^i + U 7.2 -|- 
m 3 ). If we denote (as is traditional in lunar theory) the mean angular velocity and the 
semi-major axis of this elliptic motion as n' and a\ respectively, we can write 

= Go 3 {mo + m 3 ) = Go 3 {mi +m 2 + m 3 ) . (2.25) 


For simplicity’s sake, we shall in fact consider the “Main Problem” of lunar theory in which 
the Sun is considered as moving in a circle of radius a' with the constant angular velocity 


Evaluating the derivatives with respect to X 3 in Eqs. ( 2.22 ), and using Eq. (|2.12 ), yields 


R^ = 


12- 


N.r 

n/2 


_ grav 


R 2 = m 


1 + (X 2 - Xi) 5 i2 


1 3 (N.r) 


R3 ~ mr^ (X2 - Xi) 


2+3 

5(N.r)^ -3(N.r) 

2p4 


(2.26a) 

(2.26b) 

(2.26c) 


To simplify the notation, we have written r = X 12 = xi — X 2 , N = X 3 o/r 3 o (directed toward 
the Sun), r' = r 3 o, mf^^ = Go 3 m 3 , where we recall that Gqs is the weighted average ( p.24[) , 
and 


^12 = ~ 


2 , 


(2.27) 


which agrees with our previous expressions (|1 . 1| ) or (p.3|) . In the (much smaller) term R 3 we 
have neglected the 5-modffications. 
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C. Hill’s equations of motion 


Following Euler’s second lunar theory and Hill 


0 , it is convenient to refer the 
motion of the Moon to axes rotating with the mean angular velocity of the Sun. For 
simplicity, we shall consider the main lunar problem in which (i) the Sun is considered as 
moving in a circle of radius a' with the uniform angular velocity n', (ii) the Moon moves in 
the same plane as the Sun, and (hi) one looks for a periodic motion of the Moon in the frame 
rotating with the angular velocity n!. Taking into account the lunar and solar eccentricities 
e and e', and the lunar inclination J, is expected to modify our results by terms of order 
0(e^, e'^, sin^ /) < 1%. 

With respect to the rotating frame (ex(t), ey (f)), with = N directed toward the 
Sun, the position and velocity vectors of the Moon read r = Xi 2 = Xex + Yey, v = Vi 2 = 
{X — n'Y)ex + {Y + n'X)eY (the overdot denoting d/dt). When expanding the kinetic 


terms in the reduced, relative Lagrangian (|2.21D, i.e. = |(X —nW)^ + |(F + n'X)^, 


one 


recognizes the usual Coriolis (terms linear in X and Y) and centrifugal (terms quadratic in X 
and Y) effects. The centrifugal terms can be gathered with the contribution R 2 , Eq. ( 2.26b|) , 
which is also quadratic in X and Y. The resulting Lagrangian describing the dynamics in 
the rotating frame reads 


Li 2 (X, F, X,Y) = - (x2 + (^xY - YX) + F(X, Y) . 


(2.28) 


Here, E(X, F) — Gi 2 '^o/^ T -^1 T [-^2 T \n'‘^{^X'^ + F^)] + + ■ ■ ■ — Fq + Fi + F 2 + F 3 + ■ ■ ■ 

is a time-independent potential with 


Fn = 


Fi = 


Gump 

VX^ + F 2 ’ 


n'^a' 


1 + mo/mp 


S 12 X , 


Fo = rr- 


n' 


12 


I [1 + (X2 - Xi) 512] ( 3 X 2 ) + g - (X2 - Xi) 512] (X2 + f2) I 


2 1 + mp/ms 

Fp ~ (X 2 - Xi) [5X3 _ 3 ^ ^ Y^) 

Zi Qj 


(2.29a) 

(2.29b) 

(2.29c) 

(2.29d) 


In writing out Eqs. (|2.29|) we have replaced = Gppmp by the expression = 

!{1 + mp/mp) obtained from Eq. (|2.25|) . The equations of motion corresponding to the 
Lagrangian (|2.28|) read 


,dY dF 

_2 n _ = _ 

dR dt dX ’ 

d^Y ,dX _dF 

dR ^ dt dY 


(2.30a) 

(2.30b) 


We see from Eqs. (|2.29|) that a violation of the equivalence principle has several conse¬ 
quences in the lunar theory: (i) the effective gravitational constant appearing in front of the 
Earth-Moon total mass rup = irti + m 2 , namely G 12 = G(1 + + 82 ), differs from the one 
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appearing in the theory of the Earth orbital motion, Gqs defined in Eq. (|2.24|) ; (ii) there 
is a new term, linear in X, in the Lagrangian, Fi, Eq. (|2.29b|) ; (iii) the usual tidal plus 
centrifugal potential F 2 (as well as the higher-order tidal potentials) is fractionally modified 
by ^12 7 ^ 0. The effect (i) has practically no observational consequences as, for instance, the 
“GM” of the Earth is measured much more accurately from Earth satellites (artihcial or 
natural) than from the correction it brings in the Earth-Sun interaction. The effect (ii) is 
the one discussed by Newton, Laplace and Nordtvedt, that we shall study in detail below. 
As for the effects (iii) it will be clear from the following that they are numerically negligible 
compared to the effects of Fi because the corresponding source terms in the equations of 
motion are smaller by a factor r/a' 1/400, and, moreover, the corresponding solution is 
not amplihed (as the Fi-effects) by a small divisor 1 /m ~ 12 because they correspond to 
the driving frequency 2{n — n') (instead oi n — n' for Fi). Finally, as we can also neglect 
the fractional correction (1-l-mo/ma)”^ to 812 in (p.29b|) [mo/ms ~ 1/328900], we shall keep 
Eqs. (|2.29a| ) and ( [2.29d|) and replace ( 1^.29b| ), (|2.29cl ) by 


Fi ~ , 

3 
"2 


Fn ~ . 


(2.31a) 

(2.31b) 


The contribution F 3 to the potential F (octupolar tide) generates the so-called “paral¬ 
lactic” terms in the lunar motion. Compared to the usual (quadrupolar) tidal contribution, 
they contain the small parameter r/a' ~ 1/400. Hill’s approach consists of solving hrst ex¬ 
actly (in the sense of infinite power series) the dynamics defined by the truncated potential 
Fmii = F 0 + F 2 = Gi 2 mo/r-|-|n'^X^. [The parallactic terms are obtained later by perturbing 
Hill’s main problem.] In presence of a violation of the equivalence problem we have to add 
the term Fi, Eq. ( |2.31a]) , to Fmii (“dipolar tide”!). The resulting equations of motion ( p.30|) 
read explicitly 


cPX 

dF 


2 n' 


dY 

dt 


dF 


o >dX 


Gi2mQ 

+3n'2X + n'‘^a'5i2 , 

Gi2mQ 


(2.32a) 

(2.32b) 


They admit (in the general case of a time-independent potential F) the Jacobi energy integral 


1 |'^2 ^ y2^ _ y) = G = const. 


(2.33) 


D. Iterative solution of Hill’s equations 


In spite of the apparent simplicity of Eqs. ( |2.32 ) and of the existence of the hrst integral 
( p.33| ), the corresponding dynamics contains all the richness and complexity of the three- 
body problem. Hill’s idea was hrst to hnd an exact periodic solution of Eqs. ( |2.32|) (with 
<^12 = 0). The existence (in a mathematical sense) of Hill’s periodic solution, and the 


12 































convergence of the power series (in the parameter m = n’/{n — n')) giving its explicit form, 
have been proven by Wintner |^] (see [Q for more references). The existence of such a 


one-parameter family of coplanar, periodic solutions, and the convergence of the associated 
perturbation series in m, are stable under the addition of the full series of higher-order tidal 
terms [^. We are therefore conhdent that it will still exist when one adds the “dipole tidal” 
term Fi, if 812 is small enough. 

To construct explicitly the perturbation series in m giving the periodic solutions of Eqs. 
( p.32| ) it is convenient to rewrite them in terms of new variables. Following a standard 
notation m (except that we do not introduce a separate letter for the complex conjugate 
of m) we dehne 


u = X + iY , 

(2.34a) 

u = X-iY , 

(2.34b) 

T = {n- n')t + To , 

(2.34c) 

C = e*" , 

(2.34d) 

Id d 


D = -= C— . 

i dr ^ d( 

(2.34e) 


Here, n denotes the mean sidereal orbital velocity of the periodic solution one is looking for 
(in other words the rotating frame quantities X, T, -u, u are supposed to be periodic functions 
of r with period 2tt). The parameter m = n'/ {n — n') is the remaining free parameter of the 
problem. It takes positive values for prograde orbits (going in the same sense as the Sun^: 
0 < n' < n), and negative values for retrograde orbits n < 0. With this notation the general 
equations of motion (|2.30|) read 


D^u + 2mDu 


dP 


(2.35) 


(and its complex conjugate: u u, D ^ D = —D) where F = {n — n')~‘^F = m?n'~'^F. For 
small values of m (i.e. orbits of small radius around the Farth) the periodic solution of Fq. 
( p.35| ) is of the form u ~ (const.)(C, u ~ (const.)C = (const.)(C“h It is convenient to replace 


M by a variable w which tends to zero with m. Following Liapunov and Brumberg and 


Ivanova |H, one dehnes hrst a hducial lunar semi-major axis a by writing 


GuTn-o 


(n 


n'Ya^ 


= K{m) , 


(2.36) 


where 


/ \ *^9 

K{mj = 1 -|- 2 m + -m 


(2.37) 


With this dehnition of a one introduces the variable w by 


^We shall not consider the prograde orbits with 0 < n < n' which are highly unstable; see below. 
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u = aC{l + w) , 
it = aC~^(l + w) . 


(2.38a) 

(2.38b) 


The Lagrangian C = —2a ^{n — n') ^Li 2 can be written as 

C{w, w, Dw, Dw) = DwDw + 2(m + l)wDw 

—G{w,w) + (total derivative) , 


(2.39) 


where G = 2a + {1 + 2m)(1 + w)(l + tc), and the associated equations of motion read 


dG 

D^w + 2(m + l)Dw = —7— . 

dw 


(2.40) 


The explicit form of the potential G{w,w) in our problem (i.e. when (n — n'^F = F = 
Fo + Fi + F 2 ) reads 


G{w,w) = K{m) 2{l + w) ^^‘^{1 + w) + (1 + tc) (1 + tc) +A C(l + '?^) + C ^ (1 + tc) 


.7-0-V2 


-m 


C^(l + n;)^ + r^(l + n;) 


(2.41) 


where we have introduced, instead of 5 i2, the small dimensionless parameter 


A = m^S 


^ a 

12 — 
a 


(2.42) 


The corresponding equations of motion read 


D‘^w + 2 (m + 1) Dw + K 1 + w — (1 + w) (1 + w) + AC ^ + -m^C ^ (1 + w) = 0 . 

(2.43) 

A last transformation |^], consists of separating off the square bracket multiplied by k 
in (|2.43|) its nonlinear piece, namely 

133 15 3 

Q{w, w) = (1 + w)~^^‘^ (1 + — 1 + -w + -w = -w'^ H- w'^ + -ww + O (w^, w‘^w, ww‘^, w^ 

X X rS rS 41 \ 


- 3 / 2 ' 


We can now define a /mear operator acting on {w,w), 


(2.44) 


L{w, w) = D^w + 2 (m + 1) Dw + -k (m) {w + w) , 


(2.45) 


and an effective source term (containing source terms and non-linearities) 


W{w, w) = —AC ^ {1 + w) + K (m) Q {w, w) 


(2.46) 


In terms of these dehnitions the equations of motion read 
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L{w,w) = W{w,w) , 


( 2 . 47 ) 


and its complex conjugate equation. This is the form used by Brumberg and Ivanova 
in their study of relativistic effects in the lunar motion. 

Note that the source of all equivalence-principle violation effects is the contribution 
—on the right-hand side of Eq. ( |2.46|) , with A dehned by Eq. ( |2.42|) . Even when 
restricting oneself (as we shall) to the effects linear in A, the corresponding contributions in 
the solution are quite complicated because of the interplay with the quadrupole tidal effects 
(from F 2 ), i.e., mathematically, because of the mixed term and the nonlinear 

term oc Q{w,w) in W{w,w). 

When putting back the octupole tidal effects (F 3 , Eq. (|2.29d| )) they add to the effective 
source term (|2.46|) ) the contribution 


Wr- 


3[W,W) = 


3 . 

--TT 


C(l + tc)^-F 2 C ^(1-F «;)(!-F ta)5C ^(1-Ftc) 


(2.48) 


where 


TT = {X2 


, a o m2 - mi a 

Xi)- = m -^-^ 

a m2 + mi a' 


(2.49) 


The equation ( p.47| ) can be solved by iteration: Erst, one keeps only the linear source 
terms which exist when tc = 0 , namely C~^) = hE( 0 , 0 ) = —AC“^ — [with 

the addition of 1 ^ 3 ( 0 , 0) = —| 7 r((C -|- 2 ^“^ -|- when including parallactic terms]. Second, 

one solves the linear equations L{w,w) = hE(^^((C, C”^) to get the corresponding first-order 
solution: + m‘^w ^2 + which is valid up to terms of higher order in the 

formal expansion parameters A, and n. The next step is to insert the hrst-order solution 
in the full source term W{w, w) and to collect the second-order source term C”^) 

of formal order (A -|- -|- tt)^. The corresponding second-order solution is obtained by 

solving the linear equation L{w^‘^\ = hE^^^(C, C~^)) etc.. At each stage of the iteration, 

one deals with a source term which is a linear combination (with real coefficients^) of a hnite 
number of integer powers of C and say 


w* = Wo + WiC + W_iC“^ + W 2 C" + w_2C“' + ■ • ■ + WfcC" + W.kC^ . (2.50) 


It is easy to check that there is a unique solution of the linear system L{w,w) = W*, 
L{w,w) = W* and that it is given by 

W^ = Wo + wiC + W-iC~^ + W2C^ + W-2C~‘^ H-h WkC^ + W-kC^ , (2.51) 


with 


18 


^The reality of all coefficients Wk in Eq. (|2.50 ) and Wk in Eq. ( 2.511) is easily proven by induction, 
given the reality of the coefficients in the exact definition oiW{w,w) and in the iterative solution 


Eq. (PI). 
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Wo = 


Wk = 


W-k = 


3K{m 
1 

Afc 

1 


-IVo, 

k‘^ -2{m + l)k + Wk - ^KW_k 


Al 


k'^ + 2 {m + 1) k + -k) W- 


k - -nWk 


(2.52a) 

(2.52b) 

(2.52c) 


Here Ak{= A_fc) denotes the determinant of the 2x2 system of eqnations satished by Wk 
and W-k (when k 0). Its valne is 


A, = k^ 


/c^ + 3 k — 4 (m + 1)^ 


= k^ 


k‘^ — 1 — 2m H— 
2 


(2.53) 


This determinant never vanishes {k ^ 0), bnt it takes a small valne of order m when k = ±1. 
This small divisor is one of the origins of the pecnliar amplihcation which affects both 
eqnivalence-principle-violation effects (lT*'^^(Fi) = — AC~^) and a part of the octupolar-tide 
effects (hT*'^^(F 3 ) = —|7r((C+2C“^+- ■ ■))• ^ consequence of the small divisor Ai = —2m+|m^ 
is that, when collecting from the iterative solution the contributions proportional to 
and they are found to proceed according to the powers of m^/Ai = 0 {m) instead of 
the powers of as formally expected from the structure (p.46|) p|. 

The hrst steps of the iteration can be done by hand. From Eqs. (|2.52|) the linearized 
solution (without parallactic terms) has the form 


w 


with 


) = ^ -t- w^^lc ^ , 

(2.54) 

(1) 1 2 ^ \ 

(2.55a) 

(1) 1 “1“ 2(?7i -|- 1) -|- « 

= A, ^ ' 

(2.55b) 

(1) ,9 2 

■ 

(2.55c) 

(1) 3 4 -|- 4(?7i -|- 1) -|- 2 

2 Ao 

(2.55d) 


Here, as dehned above, k = l + 2m + |m^, Ai = —2m + ^m^ and A 2 = 4(3 —2m + |m^). The 
insertion of the linearized solution Eqs. ( [^.55| ), into lF(w,w), Eq. ( |2.4(j|) , generates a 
second-order source term with the following structure: lE^^^ ~ Am^(C^^-f(C^^)-|-m^((C‘^-|-(C^^+ 
C^"^) +0{X^). Let us focus on the terms in the solution which are linear in A and contain the 
“resonant” frequencies Their source terms are found to be = ||mA (1-|-C>(m)) and 


^Following Poincare (see also Ref. |^), one can clarify the iterative process by giving a new 
name to the parameter m appearing in the second term of W, Eq. ( p.46 ), leaving unchanged the 
other occurrences of m in L and K{m). 
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w'fi = —^mX (1 + 0{m)). From Eqs. (|2.52|) the corresponding solutions can be written as 
= —f|A (1 + 0{m)), (1 + 0{m)). At this approximation we have, when 

expanding in powers of m, 


+ O 

O \ 

= - (l + Qm + O , (2.56a) 

4m ^ V // 

+ O if-i 

” 3^ +y™+ ° (™")) ■ 


E. Radial and angular perturbations due to a violation of the equivalence principle 


Let us relate the results (|2.56| ) to the radial and angular perturbations of the lunar 
motion associated with the parameter A oc 612 - The radius vector r = (X^ + and 

longitude 6 of the Moon (with respect to the rotating vector = N, i.e. with respect to 
the Sun) are such that u = re*® = a(C(l + w), where we recall that ( = e*’’. Hence 


= uu = a^(l + r(;)(l + tc) , 


^2ie 


u 


e = — = e 
u 


2iT 


1 + tc 
1 + tD 


(2.57a) 

(2.57b) 


Working linearly in A we get the following radial and longitudinal equivalence-principle- 
violation perturbations 

'1 + m}\ 

dxw 

!/ 


a 


5\6 = Am 


1 + tc 
^xw 
1 -I- tc 


(2.58a) 

(2.58b) 


At the approximation (|2.56|) we can write 5\r/a ~ (tci + tc_i) cost + {w^ + tc-a) cos3r, 
and 5x0 ~ (tci — tc-i) sin r -|- {w^ — W- 3 ) sin 3r. In the approximation, the observable synodic 
effects are entirely described by wi ± W-i. However, in higher approximations, tc± 3 ,ta± 5 , 
etc... feed down to the synodic effects in r and 6 . [Let us note in passing that, when 
averaging over time, the mean shift of the cartesian components u = X + iY = a(C(l + ta) is 
given, to all orders, by w_i alone: {X) = ata_i, (F) = 0.] 

Focussing on the contributions at the synodic frequency n — n', we get at this stage (in 
agreement with Appendix A) 



synodic 


3 A 
2 m 


14 / \ 

1 -I-m + O 


cos T , 




synodic 


= -3- 


A 


m 


1 + —m + O 


sm T 


(2.59a) 

(2.59b) 
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A straightforward, though slightly more involved, calculation allowed us to compute 
by hand the 0{m‘^) contributions to the square brackets on the right-hand sides of Eqs. 
( p.59| ). In particular, we found that the square bracket in the range perturbation, Eq. 
( p.59a|) , reads 1 -f- ym -|- -|- 0{m^). In view of the large coefficients appearing in 

this expansion, which create large corrections to the lowest-order effect (for the Moon, 
ym = 0.3773 and = 0.1445), we have decided to take the bull by the horns and to 

solve iteratively the equations of motion ( ^.47| ) to a very high order by using the dedicated 
computer manipulation programme MINIMS written by M. Moons from the University of 
Namur (Belgium) (see ^^). Some details on the application of this programme to our 
problem are given in Appendix B. Let us quote here the form of the results. Replacing A 
by its definition (|2.42|) , we see that a drops out when writing the range perturbation Sxr. 
Finally, we can write 


where 


= C (m) Si2a' cost , 
(^A0),y,odic = -C" M Aia^-sinr , 

(Jj 


C {m) 

C (m) 


3 

-m 

2 


1 + XI 


k>l 


3m ( 1 + X 

k>l 


= -m S (m) , 
= 3m S' (m) . 


(2.60a) 

(2.60b) 


(2.61a) 

(2.61b) 


The beginning of the power series S{m) entering the synodic range perturbation is 


14 

S'(m) = 1 -|- + 

6729119 


1061 


-m -|- 


2665 


48 24 

1656286531 


m -|- 


145683 

256 


-m 


2304 


-m -I- 


110592 


m® -I- 


(2.62) 


The coefficients of the series S{m) and S'{m) are given in Appendix B up to the power m^^ 
included. They are found to grow fast. The ratio between two successive coefficients Ck/ck-i 
or is found, numerically, to converge rapidly to the value 5.1254717..., thereby 

mimicking a geometric series in m/mcr with m^ = (5.1254717.. .)“^ ~ 0.195103996.... In 
the case of the Moon, with m = 0.0808489375 ... 0, this means that the series S{m) and 
S'{m) converge rather slowly, as geometric series of ratio m/mcr — 0.4144. The truncation 
to order m^^ is just enough to estimate the values of the series to the 10“® accuracy. We 
shall discuss below a method which allows us to improve this precision. We find, for the 
Moon, 


5 = 1.62201... ,5'= 1.72348 


(2.63) 


so that the full coefficients appearing in the synodic effects ( p.60|) are, respectively, C = 
0.196707..., C = 0.418025.... Finally, using the recommended value of the semi-major 
axis of the Earth orbit, a' = aEarth^ = 1.495980221 x 10^^ cm (where A denotes the 
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astronomical unit), the amplitude of the range oscillation of the Moon due to an equivalence 
principle violation is numerically found to be 

Ca %2 = 2.9427 x cm . (2.64) 


In the case where one assumes the absence of violation of the weak equivalence principle, 
i.e. (5 a = 0 in Eq. (|2.1CI|) , the result (|2.64|) gives 


((^A?^) synodic = IS.lOr^cosT meters , 


(2.65) 


if we use Eq. ( p.5|) |]33| as nominal value for the difference of gravitational binding energies. 
Our hnal result is approximately 60 % larger than the lowest-order estimate hrst derived 
by Nordtvedt in 1968 and recalled in Eqs. (p7^), ( |1.5|) and ( p78|) above. On the other 
hand, it confirms the recent finding of Nordtvedt that the interaction with the orbit’s 
tidal deformation significantly amplifies the synodic range oscillation and substantiates it 
by providing explicit and accurate expressions for the amplitude of the synodic effect. 


III. PHYSICAL DISCUSSION 


A. Resonances and instability 


We have seen in the previous section that the series in powers of m giving the amplitudes 
of synodic perturbations (|2.6U| ) appear to be close to geometric series in m/rricr with rricr — 
0.195104. This suggests the existence of pole singularities oc (/«„. — m)“^ at m = rricr — 
0.195104 in those series. Nordtvedt suggested also the presence of such a pole singularity 
at m ~ 0.2 (i.e. for a sidereal period of about 2 months) on the basis that for such a high 
orbit the driving frequency (n' in a non-rotating frame) might become equal to the rate of 
perigee advance (dzu/dt). We have substantiated and generalized this suggestion, as well as 
obtained by several independent approaches a much more precise value for rricr, namely. 


rricr = 0.1951039966... 


(3.1) 


by making use, notably, of the work of Henon on the three-body problem. To relieve 
the tedium the details of our arguments are relegated to Appendix C. Let us summarize our 
approach and our results. 

Our approach consists of putting together the (numerical) results of Henon on the 
stability of the periodic orbits in Hill’s problem, with some knowledge of the general structure 
of Hamiltonian perturbations, and a more specihc use of the analytical structure of the 
solutions of Hill’s variational equations in presence of “forcing” terms, such as the ones 
coming from the potentials Fi, Eq. ( p.29b| ), and F^, Eq. (|2.29d| ), which are neglected in 
Hill’s main problem. Our conclusions are that when m increases up to rricr, Eq. (0. 
there is a confluence of correlated singularities: on the one hand, as found by Henon, the 
free perturbations of Hill’s orbit (those not driven by any additional forces) loose their 
stability, and on the other hand, all perturbations driven by perturbing potentials of any 
odd frequency in the rotating frame (i.e. containing terms oc exp[±(2/c -|- l)ir]) develop pole 
singularities cx {rricr — rn)~^. As indicated by Nordtvedt the value m = rricr does correspond 
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to a simple 1 : 1 commensurability dw/dt = n' oi c = dl/dr = 1 in the rotating frame 
(where I = nt + e — w is the mean anomaly). Note, however, that, contrary to what happens 
in the familiar case of a harmonic oscillator, the basic frequency of the driving force does not 
need to have a 1 : 1 resonance with the natural frequency of the orbit (perigee precession) to 
generate poles oc {rricr — m)”^; the odd commensurabilities 3 : 1, 5 : 1, etc. generate similar 
poles. Therefore, both the (hypothetical) equivalence-principle-violation effects (potential 
Fi), and many of the (really existing) parallactic effects (potentials F 3 , F 5 , ...) will have 
pole singularities oc {rricr — rn)~^. Moreover, these poles are present not only in the synodic 
terms (that we concentrate upon here) but in the terms at frequencies 3 (? 7 , — n'), b{n — n'), 
.... The consequences of this situation are explored in the following subsections. 


B. Fade approximant of equivalence-principle-violation effects 


The analysis of Appendix C shows that the amplitudes of the synodic perturbations 
( p.60| ) considered as functions of m have a simple pole (but no branch point) on the positive 
real axis at m = rucr, Eq. O. and have no singularities on the negative real axis down 
to m = —1 (because of the stability of the retrograde orbits ^^) ^ This simple analytical 
behaviour suggests that the numerical validity of the power series representation ( ^.bl| ) can 
be efficiently extended by using Fade approximants, i.e. by rewriting the power series S{m), 
S'{m) truncated at order as quotients N{m)/D{m), N'{m)/D'{m) of two power series 
truncated at order m®. The explicit coefficients of the Fade approximants Q are given in 
Appendix B. We have done several checks of the conjecture that these Fade approximants 
provide a numerically accurate representation of the exact solution S{m) on the entire 
interval (—l,mcr-). First, the real zeros of smallest absolute value of the denominators 
D{m) and D'{m) are respectively found to be 0.19510399668... and 0.19510399660... in 


excellent agreement with Henon’s value (|3.1|) . Second, we found that the Fade approximants 
truncated to order md numerically agree all over the interval (—l,mcr) with those at order 
within better than 1 %. [Actually, the difference is much smaller than 10 “^ except very 
near m = —1.] We plot in Fig. 1 the Fade approximant of the coefficient C{m) of the radial 
synodic effect. 


Cpade("l) 


3 N^{m) 

2^Ds{m) ’ 


(3.2) 


over the interval {—1,rricr)- Let us note the two numerical values (using ruMoon = 
0.0808489375...; 0) 


®The value m = — 1 corresponds to very wide retrograde orbits 0 < —n <C n'. 

^It is to be noted that, thanks to the nearly geometric progression of the coefficients in many of 
the power series of Appendix B, a simpler (though less general) alternative to Fade approximants 
would be simply to factorize (1 — m/mcr)~^■ 
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Cpade(^^Moon) — 0.196707 , 
C'pade(-l) = -0.267706 . 


(3.3a) 

(3.3b) 

(3.3c) 


In Fig.l the lunar value (3.3a) is indicated by a dot. The difference between the linearized 
approximation |m (dashed line in Fig.l) and the exact value of C{m) (solid line) illustrates 
the importance of nonlinear effects in the radial synodic perturbation. 


21 



FIGURES 


FIG. 1. Coefficient C{m) of the synodic range oscillation (defined in Eq. (2.60a)) as a function 
of m = n '/(n — n'). The solid line represents the Fade approximant of C{m), while the dashed line 
represents the linearized approximation |m. The dot indicates the actual lunar value. 


C. Better orbital tests of the equivalence principle? 


Nordtvedt [^] has suggested to make use of the resonance at m = nicr to improve 
the precision of equivalence principle tests. The idea would be to put an artificial satellite 
in an orbit close to the resonant orbit (m = rricr)- From our numerical estimates, the 
resonant orbit has a sidereal period Ter = mcri^ + fTicr)~^T' (where T' = 27r/n' = 1 year), 
i.e. Ter = 1.95903 month. The corresponding “bare” semi-major axis uq = {Gmo/n?Y^^ 
is 1.68255 ao(Moon). Though interesting, this suggestion is fraught with difficulties. First, 
our results show that one must be very careful to use a sub-critical orbit (m < rrier) as 
super-critical orbits are exponentially unstable (real characteristic multiplier > 1 ). Second, 
the fact that all the parallactic perturbations (proportional to with fc > 4, i.e. 

to rn^{a/a'Y~Y develop also pole singularities at m = nier probably implies that the orbit 
becomes unstable slightly below the ideal Hill value (p.l|) P|. Moreover, the blow up of the 
parallactic perturbations amplify already large synodic effects which are known only with 
hnite accuracy. Indeed, for some years the precision of the Lunar Laser Ranging experiment 
has been limited by the accuracy with which one could theoretically compute and subtract 
the synodic parallactic perturbations 


To investigate this point we have included the 


octupole-tide perturbation F 3 , Eq. (|2.29d|) , i.e. we added the contribution IF 3 , Eqs. ( p.48|) . 


in our Hill-Brown iteration program. Our explicit results are given in Appendix B. The form 
of the radial perturbation is 


15 m 2 — mi 


mi -|- m2 


16 m 2 + mi \mi + m 2 + ms 


('^7r’")synodic 

P (m) = [K{m)]~^^'^ Q (m) 

Q (m) 


2/3 


a'P {m) cos r , 
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16 


■m + 


57599 


960 


-m + 


917401 


2880 


I 

-m -|- 


230247737 


138240 


-m -|- 


14206254151 




1658880 


The numerical value of the coefficient giving the scale of 6 .^r in Eq. (|3.4a|) is 6 .^r ~ 
28716.38P(m) COST kilometers. The Fade approximant of the series Q{m) is given in Ap¬ 
pendix B. For rough orders of magnitude estimates we can approximate (when 0 < m < mer) 
P{m) by P{m) ~ mJ^'^ {l—m/mer)~^- By comparison, the coefficient entering the Nordtvedt 
effect (|2.60a|) can be roughly approximated by C{m) ~ |m (1 — m/mer)~^- Let us dehne 


(3.4; 

(3.41 

(3.4( 


^Unpublished calculations of Henon (private communication) for a small but non-zero mass ratio 
H = mo/{mo ms) = 10“® show that the topology of the loss of stability of Hill’s prograde orbits 
is different from the ideal Hill case (/r = 0) and the same as for generic values of /r A 0. The 
difference takes place in a region of fractional size 10 “^ (~ which suggests that the actual 

is roughly 0.1 % smaller than the value (3.1). 


m. 
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the amplification factor oi the Nordtvedt effect as the ratio A{m) = C(m)/C{muoon) where 
C{m) is the coefficient in Eq. (p.60a|) . The amplification factor in the synodic parallactic 
oscillation ( |3.4a| ) will be B{m) = P{m)/P{mMoon) — {rn/mf.rY^^A{m). For an artificial 
satellite (mi -C m 2 ) one expects from Eq. (|3.4a|) that the main uncertainty in the theoret¬ 
ical value of will come from the Earth/Sun mass ratio: m 2 /m^. The current fractional 
uncertainty on this ratio is es 10“® with es — 1 |0. The corresponding uncertainty in 
is 0.07368 i?(m) centimeters. Therefore the use of a higher orbit, amplifying the Nordtvedt 
effect by a factor A{rn), will entail a correspondingly increased uncertainty on the synodic 
parallactic radial oscillation: 

5.„r ~ 0.07368 A (m) cos r cm . (3.5) 


The problem might cure itself by the fact that the ratio m 2 /m 3 will enter also the cos3r 
parallactic effects, which will be also amplified. But things might get complicated because, 
as one approaches the resonance, several frequencies become close to each other and one 
needs long data span to resolve the various frequencies and measure separately their Fourier 
coefficients. Moreover, the real motion of an artificial satellite beyond the Moon’s orbit will 
be very complex because of the combined gravitational effects of the Earth and the Moon. 
Finally, such a satellite would need to be endowed with a very high performance drag free 
system to compete with the Moon which is, naturally, drag free to high precision. 

In view of the difficulties associated with near-resonant lunar-type orbits it is worth 
thinking about other possibilities^]. Let us list some possibilities: artificial satellites around 
outer planets would be interesting in that the basic dimensionful scale factor in the synodic 
effect (|2.60a ) is a', the semi-major axis of the considered planet around the Sun. That 
would give a factor 5 for Jupiter and a factor 10 for Saturn. In either case one would need 
far enough satellites (i.e. m big enough) to have a coefficient C'(m) at least comparable 
to the lunar value (|3.3|). A second possibility is to use retrograde orbits which are always 
stable (in the Hill approximation). However, Fig. 1 shows that they give, at best, a factor 
C(—1) = —0.267706. An equivalence-principle mission consisting of a pair of artificial 
satellites around an outer planet (one prograde, one retrograde), with a laser link between 
the satellites, could improve by a significant factor upon the LLR experiment. Besides an 
improved scale factor a', the advantage of being around an outer planet is that the radiation 
pressure from the Sun is much smaller, so that the requirements on the drag-free system are 
much less stringent. 


D. Theoretical significance of orbital tests of the universality of free fall 


As we mentioned above, the Lunar Laser Ranging experiment is sensitive, through the 
synodic effect ( |2.60|) , to the sum of two physically independent contributions 


®Let us note in passing that, because of tidal dissipation, the Moon itself is, kindly, slowly receding 
toward higher orbits. However, even under the overoptimistic assumption that the present rate 
of energy dissipation continues to apply in the future, the increase in the semi-major axis of the 
Moon will be only ~ 23 % in 6 billion years (which is the expected lifetime of the Sun) ||5l||. 
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(3.6) 


<^12 — 5i — 62 'f] 


\mic^ 

= ( 5 i 2 + 4.457710-^° . 


Er^\ 


m2C^ 


The first contribution, Su = (5i — ^ 2 , is essentially equivalent to what Newton and Laplace 
had in mind when they proposed orbital tests of the universality of free fall: bodies of 
different internal compositions could fall differently. The second contribution, proportional 
to r] = 4/5 — 7 , was discovered by Nordtvedt who was assuming that the most natural 
theoretical framework in which to study possible macroscopic deviations from Einstein’s 
theory is the class of metrically-coupled theories of gravity (see e.g. 13 for a review). 
Actually, from the perspective of modern unihed theories the class of metrically-coupled 
theories of gravity seems quite ad hoc. For instance, string theory does suggest the possibility 
that there exist long-range scalar helds contributing to the interaction between macroscopic 
bodies and thereby modifying the standard predictions of general relativity. However, all the 
scalar helds present in string theory have composition-dependent couplings for very basic 
reasons (for a discussion of general theoretical alternatives to Einstein’s theory and the types 


of composition-dependent couplings they might exhibit see |]35[| ). 

Recently, a mechanism has been proposed by which some of the scalar helds of string 
theory might survive in the macroscopic world as very weakly coupled long range helds 
(see also |^). In the model of Ref. [M the surviving scalar held(s) modify the ob¬ 


servational consequences of general relativity in several ways: (i) they violate the “weak 
equivalence principle” {Sa 7 ^ 0 ) because of the composition-dependence of their couplings to 
matter; (ii) they modify the post-Newtonian {0{l/c^)) ehects in essentially the way which 
is parametrized by the Eddington parameters^ (3 and 7 ; and (hi) they induce a slow time 
variability of all the coupling constants of Nature: G, a, aweak, ■ ■ ■ ■ The point we want 
to emphasize here, because we think it is generic, is that all those modihcations of general 
relativity are related, because they derive ultimately from the couplings of the same held. 
In particular, the hrst term on the right-hand side of Eq. (|3.6|) is related to the second. 

More precisely, in the model of Ref. [^, we have, for an individual atom labelled by A, 
the link 


h ^ -0.943 X 10“®7(E/M)a (3.7) 

where E = Z{Z — 3)/{N + ZY^^ is associated to the electrostatic interaction energy of the 
nucleus of the atom, and where M denotes the mass of A in atomic mass units. We believe 
that the structure of this link between 5a and ^{E/ M)j^ is generic in string-derived models, 
because it follows from a basic physical feature of the massless scalar helds (“moduli”) 
present in string theory, namely that they dehne the values of the gauge coupling constants. 


^^This comes from a feature of their couplings which is deeply rooted into the structure of QCD 
and the consequences it has for the mass of atoms; see pages 550-553 of |@. 

We neglect here contributions proportional to the ratios (baryon number)/(mass) and (neutron 
excess)/(mass) which tend to be subdominant, even for moderate Z differences |36||. 
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Even the magnitude of the numerical coefficient shoud be somewhat generic. Indeed, its 
analytical expression —lasaAc/A^g (in the notation of [^) shows that it is determined by 
basic physical facts or assumptions: fractional smallness of electrostatic nuclear contributions 
( a^a ~ 0.770 X10“^), unihcation of gauge coupling constants (A^ — 1) and QCD conhnement 
(A„3 ~ ln(Astring/(a.m.u.)) ~ 40.8). 

We have also the model-dependent link (3 ~ —10.2^7, where the dimensionless theory 
parameter n is expected to be of order unity, [k > 0 denotes here the curvature of a coupling 
function around a minimum and should not be confused with the notation K{m) used above.] 
These links indicate that, in the Earth-Moon case, the gravitational binding contribution 
to 6 i 2 is numerically negligible compared to the composition-dependent term Si 2 - Indeed, 
using Eq. ( p.7|) and the compositional difference between the Earth and the Moon (i.e. the 
difference between an Earth iron core of mass 0.32m2 and a silica-dominated Moon [^]), we 
hnd 5 i 2 — 0.32((5si — (^Fe) — 2.75 x 10“®7, while the gravitational binding energy contribution 
is 4.457 X 10-1° = -4.45(40.8k + 1) x 10' ■i°7. From the point of view advocated here, the 
conclusion is that LLR data give us a very precise test of the weak equivalence principle. 
The loss of a Nordtvedt-type direct test of the combination 7 = 4/? — 7 is compensated by 
the theoretically expected link 812 — 812 — 2.75 x 10-°7 which gives an extremely good limit 
on the effective Eddington parameter of the considered scalar model. More precisely, the 
observational limit[^ 812 = (—3.2 ± 4.6) x 10-i° recently derived from LLR data 
translates into the following observational constraint on 7: 


50 , P 


7 = (-1.2 ± 1.7) X 10 


-7 


(3,8) 


The recent laboratory tests of the weak equivalence principle give comparable results. Using 
the experimental limit ^BeCu = (—1.9 ± 2.5) x 10“^^ [Q and the theoretical formula ( p.7|) 
(which yields (fBeCu = 2.41 x lO-^y) we hnd 


7= (-0.8 ±1.0) X 10 


-7 


(3.9) 


Impressive as these limits may seem. Ref. gives a motivation for pushing equivalence 


principle tests further because this reference estimates that a natural range for the coupling 
parameter 7 is 10“^° ^ —7 ^ 10“^°. [Note, however, that if the theory parameter k is 
of order 1/40 (which corresponds, in the notation of [^|, to k ~ 1) larger values of —7, 
of order 10“^, are expected, in agreement with |^.] In this connection, let us mention 


that the LLR CERGA team plans to improve the precision of the ranging down to the 2-3 
millimeter level for normal points [C. Veillet, private communication]. Extracting 812 at this 
level will necessitate to improve the modelling of the solar radiation pressure effects which 


are currently believed to contribute a synodic range oscillations of approximately 0.3 cm [25 


^^This was obtained from partial derivatives of the numerically integrated equations of motion, 
and therefore independently of theoretical estimates of the value of the coefficient C{m) in Eq. 
(|2.60aD . 

^°We do not take into account here the theoretical constraint that 7 < 0 in all scalar models. 
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If this can be done, the LLR experiment will reach the level J12 ~ 5 x 10“^“^ corresponding to 
the level 7 ~ 10“®. It seems that significant progress in testing the equivalence principle will 
require space missions: either a low Earth orbiting artihcial satellite dedicated to testing 
the weak equivalence principle, as the STEP (Satellite Test of the Equivalence Principle) 
mission, or, possibly, some type of orbital test such as the one suggested in Ref. and the 
ones suggested above. 
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APPENDIX A: TRADITIONAL LUNAR PERTURBATION THEORY 


As a check on the lowest orders of the Hill-Brown calculations presented in the text, we 
have also investigated the mixing between equivalence-principle-violation effects and tidal 
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effects by means of the traditional lunar perturbation theory of de Pontecoulant 
The equations of motion corresponding to the Lagrangian L = + n/r + R read 




dY r _dR 
dR ^ dr 


(Al) 


Here, r = X12 = xi — X2 is the position vector of the Moon with respect to the Earth (in 
an inertial, non-rotating, coordinate system), /i = Gi 2 {mi + m 2 ) and R = Ri + R 2 + R 3 + 
■ ■ ■ is the total potential perturbing the Keplerian motion of the Moon around the Earth. 
[This corresponds to Eqs. (|2.21D , (p.22|) .] We consider the coplanar problem for which it is 
enough to solve for the radius r = |r| and the longitude v (polar angle). Decomposing the 
acceleration into radial and longitudinal components leads to 


r — = —fir ^ -|- dR/dr , (A2a) 

d{r^v)/dt = dR/dv . (A2b) 


de Pontecoulant’s method uses the longitudinal equation ( |A2bD but replaces the radial one 
( |A2a|) by the “virial” equation dealing with the second time derivative of r^. The basic 


equations are then written as 


1 d^ 
2 dR 



p ^ 

r ttc 
he 


P . 

Lfrf — 


(A3a) 

(A3b) 
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where Oc and he are some integration constants and where the transformed source term in 
the radial equation is 



(A4) 


In the hrst form of P, {d/dt)i denotes a time derivative taking into account only the variabil¬ 
ity due to the time-dependence of the coordinates of the Moon: {d/dt)i = rd/dr -|- vd/dv. 
The second form of P is obtained by taking into account the time-dependence of R upon the 
Sun’s coordinates, and assumes that the Sun moves on a circular orbit (r' = 0, h' = n'). It 
is very useful to notice that if P = Rp where each contribution Pp(r, v — v') has a radial 
dependence oc r^, then P = Pp with 


Pp — (p -|- 2) Rp 2ti J 


dt 


dRp 

dv 


(AS) 


Note that the use of the suffix p is consistent with the notation Pi, P 25 R 3 of Eqs. ( 2.22 ). 

First-order perturbation theory is very easy. Let us consider a general term R^^ = 
Ar^ cos q{v — v'), perturbing the zeroth-order (circular) solution r = vq = nt + e (with 
zeroth-order integration constants Oc = a, he = no?, and the link = p). Inserting 
the perturbed solution r = a + 5r, v = Vq + 5v into Eqs. (^) (with perturbed integration 
constants ae = a + 5a, he = h + 5h) yields 


{5r) + n^5r = pAaF ^ cos [q {n — n') [t — p)] + n^5a , 

(i n \ 

— {5v) -|- 2—5r = - - -AaF~‘^ cos [q in — n') it — p)] + a~^6h , 

dt a [n — n') 


(A 6 a) 

(A 6 b) 


where 

n' 

p = p + 2 + 2 -= p -|- 2 -|- 2 m . 

n — n' 


(A7) 


Here, as in the text, we shall use m = n'/{n — n') as the small expansion parameter of 
perturbation theory. The solution of Eqs. ([A^ is (with r = vq — Vq = (n — n'){t — to)) 


6 r 

6 v 


P 

iP — q^{n — 
1 

q{n — n') 


— r 7 :AaF ^ cos qr + 5a , 
n')^ 

2 pn 1 


— q^{n — n')"^ n — n' 


AaF ^ sin qr -|- {5n)t -|- 5e , 


(A 8 a) 

(A 8 b) 


where 5n = —2na ^5a -|- a ^5h. Eqs. (|AB|) 


have been written assuming q ^ 0. They 


take a different form when q = 0. It is traditional to keep n hxed throughout the 


approximation process (and therefore numerically equal to the observed mean motion), and 
to dehne a by tPa^ = p. Then 5h is computed in terms of 5a so that 5n = 0. Finally, one 
must make use of the original radial equation ([A2a|) to determine 5a. 
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We are especially interested in the case where the perturbing potential R is the sum of 


the quadrupole tide -R 2 , Eq. (|2.26b|) , and of a term with frequency g = 1: 


R — R 9 


TD /2 2 

Kp = n r 


13 

- + - cos 2 (n 
4 4 


+ Ar^ cos {v — v') . 


(A9) 


In the case of the equivalence-principle-violating term (|2.26a|) the perturbing term in 
cos{v — v') has p = 1, while the octupolar tide (|2.26c|) , for which p = 3, contains a perturbing 
term in cos(n—v') (that we focus on) and a term in cos 3(n—u'). At the fractional order 0{m‘^) 
beyond the hrst-order solution (|A8|) , the cos3(n — n') term mixes with the quadrupole tides 
oc cos2(n —n') to generate the frequency g = 1. As we work here only at the fractional order 
0 {m) beyond the hrst-order solution, we do not need to study the effect of the cos3(r; — v') 
term. The hrst-order solution corresponding to Eq. (^) reads 


r = a -|- 6^r + 6‘^r -|- S^r , 
V = nt + e + 6^v - 1 - 6^v , 


(AlOa) 

(AlOb) 


where the superscripts on 6 indicate the value of the frequency g. To sufficient accuracy for 
our purpose we have from Eqs. (|A8| ) 

(Alla) 
(Allb) 
(Allc) 

(Alld) 

yit — If)- 

with 


5^r -\- 6 ‘^r = —rn^a -|- cos 2r^ -|- 0 {w?) , 

5‘^v = — m? sin 2r -|- 0 {m^) , 
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first-order 


= a 


_ ^first—order 
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(n — n')2 
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(n — n'y 


aF ^ cos r 


aF ^ sin r , 
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'first—order 


p -|- 2 -|- 2 m p + 2 


2 m + m^ 


2 m 


1 + 


2 — p 
2(2+ p) 


m 


+ O [mF'^ 


a 


first-order 


= -2 


(p -I- 2 -I- 2m) (1 -I- m) 
2 m -|- m^ 


+ 1 = -2 


p + 2 

2 m 


1 + 


4:+P 
2(2+p) 


m 


+ O (m'^^ 


(Al2a) 

(Al2b) 


Note that the constant term —mFa/Q in Eq. ( |Alla|) depends on 5a = Oc — a and must be 
determined by having recourse to Eq. (|A2a|) . Note also that the small denominator present 
when g = 1, rF — {n — n'Y = (n — ? 7,')^[(1 -|-m)^ — 1] = (n — ? 7 ,')^[ 2 m-|-m^], causes the synodic 
ehects to be of order 0{m~^A) instead of the usual order 0{A) valid when g 7 ^ 1. One of the 
ehects of this small denominator is to have a simple, approximate link between the radial 
and longitudinal synodic oscillations: (^fi^st-order _ _ 2 (^first-order|-]^ _|_ ( 9 (^ 177 ,)]. Another ehect 
is that only the leading terms in Eqs. (A12|) are correct. The 0{m) fractional corrections 
are modihed by higher iterations as we are going to see. 

When proceeding to the next iteration, several effects must be taken into account. On 
the left-hand side of Eq. ( |A3a|) one must keep the terms of order (5r)^, while on the right- 
hand side one must include the change 5P of the source term P{r,v — v') induced by the 
hrst-order solution (|A11|). This leads to the following equation for 6 r = r — a 































= ^(^ 0 , Vo - v'o) + 6PeS , 


(Al3a) 

(Al3b) 


When computing the synodic effects with fractional accuracy l + 0{m), the computation 
of 6 P is simplified by several circumstances. Because of the amplification 0{m~^A) of the 
first-order synodic effects, one finds that it is enough to compute in Eq. (|A5|) the change of 
the first contribution, (p -|- 2)i?p, for p = 2 and under the synodic variations 6 ^r, This 
yields 


m 


d 


d 


synodic 


= 4 (<5i?2)sy„odic = \ (+ 3 cos 2 (n - n')) 


= 11 n'^a^Scosr [1 -I- C> (m)] , 


dv 


synodic 

(AM) 


where B = {2m)~^{p + 2){n — n')~‘^Aa^~‘^ denotes the leading value of the fractional synodic 
range oscillation {6^r/a ~ Boost, ~ —2i?sinr). We need also to extract the synodic 
piece of (hr)^ = (S^r P S^r -|- S^r)^ coming from the mixing between -|- S^r ~ —m^a(| -|- 
cos2r) and ~ aBcosr: (^?")synodic — —§a^v 2 ^B cost[ 1 P 0{m)]. Finally, the synodic 
piece of the second-order effective source term for de Pontecoulant’s radial equations (|A13|) 
is obtained as 


s'p,„ = + 


P 


(n — n 






synodic 


= 9 n'^a^B cos t[1 P O (m)] . 


(A16) 


The corresponding solution reads 



second—order 


~ 9 


n 


/2 


m? — {n — n'Y 


9 

~ -mB COST . 
2 


B COST 


(A16) 


When turning to the longitude equation (A3b|) one finds also some simplifications: the 
change of the source term oc J dt dR/dv is of order x m~^A = 0{mA). The corresponding 
term in the solution is not amplihed by a small denominator and is therefore negligible 
compared to the precision m x m~^A = 0{A) we are aiming for. It is therefore sufficient 
to integrate the equation d/dt ~ —The hnal result 

can be very simply expressed in saying that the second iteration leads to multiplying the 
hrst-order synodic perturbations ( |AllcD , (|AlldD by the common factor 1 -|- |m -|- 

In conclusion, the mixing between the quadrupole tide R 2 and some synodic-frequency 
perturbation potential (in which we factorize an effective gravitational mass of the Sun, 
Gm' = 


(A 


P/synodic 



(A17) 


29 



















leads, when neglecting non-linearities in the dimensionless parameter (3 (which should not 
be confused with its post-Newtonian homonym), to the following synodic oscillations 




= Cr(3rn^ 


p-2 


COST 


6^v = C^l3m? 


p-2 


where 


a = 


9 , . 

1 -I- -m + O [rn?j 


10 -|- 4p 

2 + p 
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smr 


first-order 


(Al8a) 

(Al8b) 


2 + p 

2 m 


1 + 


-m 


+ O [m?^ 


C,= 
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,2 + p 

2 m 


1 -|- -m + O (m'^'^ 
= -2 
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first-order 
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1 + ^ + O (m^) 


2 + p 


(Al9a) 


(Al9b) 


The two cases of interest are: (i) an hypothetical violation of the equivalence principle 
in which (comparing Eq. (|A17|) with Eq. (|2.26a|) ) 


P — 1 ) /3eP — <^12 , 

and, (ii) the octupolar tide (“parallactic effects”), Eq. (|2.26d) , with 

1 Gtti^ / r \ ^ 

R,c^-{X2-X,] 


whose synodic piece has 


(A20) 


a’ \a' 

X [3 cos (n — n') -|- 5 cos 3 (v — n')] , 


p — 3 , Ppar — - {X 2 — Xi) 


(A21) 


(A22) 


The results (|A18|) - (|A22|) agree with the (much more accurate) Hill-Brown-type results 
given in the text and in Appendix B. 


APPENDIX B: HILL’S EQUATIONS MORE ON THE ITERATIVE SCHEME 
OF SOLUTION AND NUMERICAL RESULTS 


In this Appendix, we explain in detail the iterative scheme we employed for solving 
Hill’s equations ( |2.47|) with the source terms { \2A^ ), and also with the parallactic perturba¬ 
tion (|2.48|) . We also give tables of the obtained solution for several physically interesting 
quantities. Obviously, one can envisage several iterative methods for solving the considered 
equations. We do not claim that the scheme we adopted is the optimal one, but we found 
it suitable from the point of view of memory and computing time requests. Thanks to our 
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use of the dedicated algebraic manipulator MINIMS we could obtain the solution to a very- 
high order in the formally small parameters. In what follows, we shall present the solution 
for the perturbation of Hill’s variational orbit related to the equivalence principle violation. 
Exactly the same scheme applies in the case of the parallactic perturbations. 

Keeping the notation of Sec. II.D, notably L{w, w) for the linear operator (p.46| ), we have 
to solve 

L{w,w) 

3 

— -m‘^'w(~‘^ + K{m)Q{w,w) , (Bl) 

where 

Q {w, w) = (1 + (1 + - 1 + . (B2) 

The nonlinear function Q can be written as 

Q (w, w) =J2xk (w^ + {‘2k + 1) w*’) 

k>2 

+ ^ {2k + l)xjXkW^w^ . (B3) 

j,k>l 


Here, Xk = 


- f -^/2 


k 


and {2k + l)xk = 


-3/2 

k 


are binomial coefficients. 


As in Sec. H.D we look for a formal power series solution of Eq. 

W = tyB) + -u;® y^(3) -I- . . . 


(B4) 


and similarly for the complex conjugate. The superscripts on the consecutive terms in ( [B4|) 
refer to corresponding orders in the combined formal small parameter (A + m?). Keeping 
track of the orders in this formal small parameter, we decompose the nonlinear source Q as 

Q(u;,u;) = g(2)+g(3) + ... ^ 


where the individual terms include symbolically 

gd) ^ ^ (coefficient) . (B6) 

ja-\-kb=i 

For any particular value of i in (|B(j|) , Qh) jg given by a finite number of terms which depend 
only on the knowledge of for a < i. Although the procedure of breaking Q{w, w) into a 

sum of equal-order terms might seem laborious, it is relatively easy to be programmed 

using a well suited algebraic manipulator such as MINIMS. One can introduce a formal index 
which conserves the order of a particular term and manipulate it like any other variable. 

The heart of our iteration scheme consists of the following infinite system of differential 
equations 
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l(w^’^\w^^^) 


-k-^ - Irn^C^ , 

— ^m‘^w^'^\~'^ + (w^^\ , 

. . . , 

-|mV‘-')C“" + kQW (w<'>, w'^K .... u.l'=-‘>) . 


(B7) 


It is easy to verify, that the generic form of the terms in the right-hand sides of Eqs. (J 
reads H4C^ + ^-kC~^ presented in (|2.50|) . The nniqne inversion of the linear operator 
L on the left-hand sides of Eqs. (0) is given by formulas ( p.51| )-( |2.52|) . Notice also that 
suppressing the equivalence-principle-violation term (A = 0 in (|B7|) ), one recovers a system 


of equations for constructing the usual variational periodic orbit. 

In the preceding scheme, we consider A and m? as two comparable “small” parameters. 
However, in practice, the order of the A parameter associated with the studied violation of 
the equivalence principle is numerically much smaller than (which can be as large as 
1). As a result, we restrict the generality of our solution by keeping only the first order in 
the parameter A. This truncation allows a clear separation in the interpretation of the odd- 
and even-power terms in the ( variable of the final solution for w (and w): (i) the even- 
power terms (oc never contain the perturbation parameter A and fully reconstruct Hill’s 
variational solution, (ii) the odd-power terms (oc s-re all of the hrst-order in A (but 

they are coupled to the “background” variational terms through an infinite series of powers 
of m?). We thus simultaneously obtain Hill’s variational solution and its A-perturbation by 
hltering the various powers of (. This is a particularly important circumstance, because the 
series giving the variational solution enters the definition of several studied quantities such 
as the radial or longitudinal perturbations of the lunar orbit by the equivalence-principle- 
violation terms (see Eqs. (p.58|) ). 

Once the iterative scheme is set up and the numerical program debugged, we can obtain 
the solution of our problem up to an arbitrary order. The limits of the solution are then 
given mainly by the computer power. A minor limit comes from the fact that the MINIMS 
algebraic manipulator works with double precision real coefficients (16 digits) [^. During 
the manipulation of the series, one thereby accumulates round-off errors. However, we have 
checked that this restriction is not signihcant for our workQ. 

In the rest of this Appendix we shall present tables of the numerical coefficients achieved 
by the previous algorithm for different series introduced in the main text of the paper and 
related to physical quantities. 

Let us start with our solution for the variational curve (A = 0). Tables I and H give the 


We have performed a lower order solution in using the modified version of the distributed 
MINIMS manipulator which accepts the quadruple precision (32 digits) for the coefficients of the 
series and compared it with the double precision one. 
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coefficients Wij of the double series expansion of w:w = where j = 0, ±2, ±4,.. 


fc = 2,3,4,... and |j| < k {j labels the rows and k the columns). Contrary to the method 
of Ref. 1^ the ^-independent term is not fixed to unity. However, due to the choice of the 


fiducial semimajor axis a, dehned in Eq. (p.36|) , it starts only at the power m^. 
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TABLES 


TABLE I. Coefficients Wjk of the double series giving Hill’s variational curve: w = J2 


j labels the rows and k the columns. 



2 

3 

4 

5 

6 

-6 

- 

- 

- 

- 

0.005208333333333333 

-4 

- 

- 

0.0000000000000000 

0.03593750000000000 

0.1245833333333333 

-2 

- 1.187500000000000 

- 1.666666666666667 

- 1.194444444444444 

- 0.5185185185185185 

- 0.8265365788966050 

0 

- 

- 

0.6210937500000000 

- 0.6770833333333333 

- 0.6892361111111111 

2 

0.1875000000000000 

0.5000000000000000 

0.5833333333333333 

0.3055555555555556 

- 0.1615849247685185 

4 

- 

- 

0.0976562500000000 

0.4182291666666667 

0.8484722222222222 

6 

- 

- 

- 

- 

0.06778971354166667 




TABLE II. 

Continuation of Table I. 



7 

8 

9 

10 

11 

-10 

- 

- 

- 

0.002728271484375000 

0.02314830588154192 

-8 

- 

0.003743489583333333 

0.02753602458235899 

0.09640929637552359 

0.2138120808825399 

-6 

0.03477027529761905 

0.1040162627551020 

0.1880280591858106 

0.2403927412241086 

0.2416131234195372 

-4 

0.1956215277777778 

0.2153863326461227 

0.2955465183221726 

0.3769043026659738 

- 0.1199168976723571 

-2 

1.205252137988683 

7.122740269204389 

15.70605739258259 

19.78110897697562 

14.66018287215882 

0 

- 1.605902777777778 

- 0.5700574333285108 

0.3975170245386445 

1.670582431118184 

7.457045029767965 

2 

- 1.034723427854938 

- 3.070661490483539 

- 6.516200630679869 

- 9.530344734021654 

- 7.910677339637028 

4 

1.038887731481481 

0.7041097601996523 

- 0.5710030937052194 

- 4.063905605879704 

- 11.57027457856626 

6 

0.3898297991071429 

1.087330552012472 

1.926893073875549 

2.301395728392179 

1.149297956477343 

8 

- 

0.05397033691406250 

0.3869401996638499 

1.360005033209030 

3.095895579274965 

10 

- 

- 

- 

0.04656556447347005 

0.3988932608444074 
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We then give the lunar orbit perturbations due to a hypothetical equivalence principle 
violation (terms linear in A). 

Tab. gives the coefficients of the series in powers of m giving tc_i, i.e. the coefficient 
of in the Laurent expansion of S\w{Q, after factorization of the leading term |(A/m) 
(see Eq. ( |2.56b| )). The second column ~ pk - gives the numerical value of (in %) for 
the lunar orbit (m = rriMoon = 0.0808489375...). The last column - - gives the ratio 

Ck-i/ck of the successive coefficients of the series (the same structure is conserved also for 
Tables and |VlTHX|). 


TABLE III. Coefficients of in the lunar orbit perturbation rc-i. The mark < signifies that 
the value is smaller than 0.001 %. 


k 

Cfc 

Pk 

Tk 

0 

1.000000000000000 

- 

- 

1 

5.111111111111111 

41.323 

0.1956521739 

2 

25.52777777777778 

16.686 

0.2002176279 

3 

129.5277777777778 

6.845 

0.1970834227 

4 

663.1076388888889 

2.833 

0.1953344679 

5 

3400.509837962963 

1.175 

0.1950024174 

6 

17434.56978202160 

0.487 

0.1950440923 

7 

89366.97811374742 

0.202 

0.1950896198 

8 

458049.7173103071 

0.084 

0.1951032273 

9 

2347711.432406117 

0.035 

0.1951047778 

10 

12033105.92287766 

0.014 

0.1951043602 

11 

61675313.45886662 

0.006 

0.1951040902 

12 

316115045.4115870 

0.002 

0.1951040115 

13 

1620238694.871019 

0.001 

0.1951039970 

14 

8304487501.579596 

< 

0.1951039958 

15 

42564415192.16274 

< 

0.1951039962 

16 

218162702733.8958 

< 

0.1951039965 

17 

1118186744062.939 

< 

0.1951039966 
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Tab. gives the coefficients Ck of the series in powers of m giving wi, i.e. the coefficient 
of C in the Laurent expansion of S\w, after factorization of the leading term — |(A/m) (see 
Eq. (p.56a|) ). 


TABLE IV. Coefficients of in the lunar orbit perturbation wi. Notation as in Table III. 


k 


Pk 

rk 

0 

1.000000000000000 

- 

- 

1 

6.000000000000000 

48.509 

0.1666666667 

2 

29.62500000000000 

19.365 

0.2025316456 

3 

147.5000000000000 

7.795 

0.2008474576 

4 

749.9427083333333 

3.204 

0.1966816910 

5 

3843.245659722222 

1.328 

0.1951326495 

6 

19711.47459129051 

0.551 

0.1949750457 

7 

101057.1821729118 

0.228 

0.1950526837 

8 

517989.9317152261 

0.095 

0.1950948773 

9 

2654938.425605692 

0.039 

0.1951043108 

10 

13607759.81112155 

0.016 

0.1951047389 

11 

69746091.35545381 

0.007 

0.1951042639 

12 

357481499.9590677 

0.003 

0.1951040581 

13 

1832261194.517921 

0.001 

0.1951040065 

14 

9391202741.857013 

< 

0.1951039973 

15 

48134343317.40093 

< 

0.1951039963 

16 

246711211454.4727 

< 

0.1951039964 

17 

1264511316090.325 

< 

0.1951039966 
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Some general properties of the presented results, already commented upon in the main 
text of the paper and in Appendix C, are: (i) a fast (“geometric-like”) growth of the coef¬ 
ficients Cfc of the series with a surprisingly rapid approach of the rfe(= Ck-i/ck) ratio to the 
nicr value of about 0.1951039966..., (ii) a very substantial contribution of the non-linear 
terms {k > 1) in the series for the lunar value of Hill’s parameter rriMoon- 


TABLE V. Coefficients of the S{m) series yielding the radial perturbation, with synodic (“r”) 
period, of the variational curve due to an equivalence principle violation. 


k 


Pk 

rk 

0 

1.000000000000000 

- 

- 

1 

4.666666666666667 

37.730 

0.21428571429 

2 

22.10416666666667 

14.449 

0.21112158341 

3 

111.0416666666667 

5.868 

0.19906191370 

4 

569.0742187500000 

2.431 

0.19512686221 

5 

2920.624565972222 

1.001 

0.19484675483 

6 

14976.54921694155 

0.418 

0.19501318519 

7 

76767.66017493731 

0.173 

0.19508930170 

8 

393469.7706768071 

0.072 

0.19510434065 

9 

2016707.919972300 

0.030 

0.19510498609 

10 

10336561.13767503 

0.012 

0.19510433819 

11 

52979731.33709500 

0.005 

0.19510406861 

12 

271546096.6185635 

0.002 

0.19510400627 

13 

1391801807.779853 

0.001 

0.19510399764 

14 

7133640668.991231 

< 

0.19510399701 

15 

36563272849.28028 

< 

0.19510399680 

16 

187404017729.6833 

< 

0.19510399666 

17 

960533976540.2227 

< 

0.19510399664 
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TABLE VI. Coefficients of the S'{m) series yielding the longitudinal perturbation, with syn¬ 
odic (“r”) period, of the variational curve due to an equivalence principle violation. 


k 


Pk 

rk 

0 

1.000000000000000 

~ 

- 

1 

5.333333333333333 

43.119 

0.18750000000 

2 

26.39583333333333 

17.254 

0.20205209155 

3 

132.9166666666667 

7.024 

0.19858934169 

4 

678.3172743055553 

2.898 

0.19595058493 

5 

3478.588686342593 

1.202 

0.19499783834 

6 

17838.85192117573 

0.498 

0.19500070418 

7 

91447.83453394253 

0.206 

0.19507134326 

8 

468721.9555645895 

0.086 

0.19510038616 

9 

2402405.984973695 

0.035 

0.19510522305 

10 

12313419.50072995 

0.015 

0.19510469735 

11 

63112018.54675127 

0.006 

0.19510419385 

12 

323478806.4318246 

0.003 

0.19510402936 

13 

1657981435.159388 

0.001 

0.19510399789 

14 

8497936880.202027 

< 

0.19510399507 

15 

43555934597.05303 

< 

0.19510399579 

16 

223244707500.3778 

< 

0.19510399635 

17 

1144234415604.679 

< 

0.19510399657 
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Tables HH give the coefficients of the power series S{m) and S'{m) related, respectively, 
to the radial and longitudinal perturbations (with synodic periodicity) of the lunar orbit due 
to a hypothetical violation of the equivalence principle. The coefficients of the Fade approx- 
iniants of those series are given in Tab. More precisely, we denote Ns{m) = 
and D^{m) = the Fade approximants for the S{m) series (see Eq. (|3.2|)), and 

N'^{m) = 3.nd D'^{m) = the Fade approximants for the S'{m) series. The 

denominator polynomials Ds{m) and D'^{m) have 0.1951039967... and 0.1951039966..., 
respectively, as real roots. 


TABLE VII. Coefficients of the Fade approximants of order 8 of the radial and longitudinal 
perturbation series, with synodic (“r”) period, of the variational curve due to an equivalence 
principle violation. 


k 


bk 


b'k 

0 

1.000000000000000 

1.000000000000000 

1.000000000000000 

1.000000000000000 

1 

-4.679659654425508 

-9.346326321092175 

-2.948114664716526 

-8.281447998049859 

2 

9.706375765630793 

31.21839859739427 

2.472472675594219 

20.24436199852680 

3 

-12.09649460259615 

-62.23126666796116 

-0.035120215085467699 

-22.32599642537073 

4 

9.264485113438459 

78.37781038863183 

-0.8074761600241968 

6.322888291078106 

5 

-0.8206525542778040 

-59.42770492795888 

1.129701196880901 

4.768342316076093 

6 

-4.562346877625465 

5.529574170471975 

1.916287015273362 

13.88318783117331 

7 

0.6450458113816286 

30.28114298289518 

-3.678767273115525 

-38.01361281338056 

8 

0.9269455883562498 

-30.33613801897244 

-3.906968891383218 

34.97077854522613 
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Tab. |VIII| gives the coefficients of the S"{m) series yielding the equivalence-principle- 
violation perturbation of the lunar orbit with one-third of the synodic period (“3r”). This 
series is dehned by 


= ^"(mjiija'cosST , 


(B8) 


with 


C" (m) = —m^ 
^ ’ 32 


1 + E 


Cfjn 


k>l 


k 1 = gn 

32 


m) 


(B9) 


Note that the third-synodic effect ( |B8|) is 0{m?) smaller than the synodic effect ( p.bOl) . 
For instance, if we assume ^12 = 0 the result (|B^) gives numerically (5A?")third-synodic — 
6.627] cos 3r centimeters for the lunar orbit. Although there is probably no practical use of 
this higher frequency excitation (and of the others with frequencies 5r, 7r, ...), two points 
are to be mentioned: (i) a signihcantly smaller amplitude of the effect, (ii) the persistence 
of the pole singularity aX. m = rricr for these odd-multiples of the basic synodic frequency as 
discussed in Appendix C. 


TABLE VIII. Coefficients of the S"{m) series yielding the radial perturbation, with 
third-synodic (“3r”) period, of the variational curve due to an equivalence principle violation. 


k 

r" 

Pk 

rk 

1 

1.000000000000000 

- 

- 

2 

6.450980392156863 

52.155 

0.15501519757 

3 

32.60457516339869 

21.312 

0.19785506665 

4 

163.0646514161220 

8.618 

0.19994876192 

5 

828.6388478122731 

3.540 

0.19678615340 

6 

4240.758469846889 

1.465 

0.19539873674 

7 

21736.94075865620 

0.607 

0.19509454053 

8 

111428.4766430605 

0.252 

0.19507527531 

9 

571159.5105168029 

0.104 

0.19509169434 

10 

2927510.871437882 

0.043 

0.19510073083 

11 

15004916.03283528 

0.018 

0.19510344910 

12 

76907280.99332377 

0.007 

0.19510397246 

13 

394186033.1895487 

0.003 

0.19510402327 

14 

2020389145.708269 

0.001 

0.19510401450 

15 

10355446785.91400 

< 

0.19510400541 

16 

53076547791.36041 

< 

0.19510400011 

17 

272042338525.2278 

< 

0.19510399771 
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Tab. |I^ lists the coefficients of the Q{m) series giving the radial parallactic inequality of 
the lunar motion {Q{rn) = 1 + J2k>i as dehned in Eqs. ( p.4|) of the text. Similarly to 

the treatment of the equivalence-principle-violation lunar perturbation we improved on our 
solution by using Fade approximants. Tab. ^ yields the coefficients of the corresponding 
polynomials. We also computed the corresponding lunar parallactic inequality in longitude. 
As a partial check on our results we have compared the latter with the result by Deprit, 
Henrard and Rom |^. When substituting the current recommended values of the mass 
constants and the Earth semi-major axis we obtain 125.” 438 for the amplitude of the par¬ 
allactic inequality in longitude when truncating our series to the power m7. This value is to 
be compared with 125.” 4201 reported in Ref. [^. We believe that the origin of the minor 


discrepancy between those results lays in the slightly different values of the astronomical 
constants employed by Deprit, Henrard and Rom at the end of the sixties. 


TABLE IX. Radial perturbation of the variational curve due to parallactic terms with synodic 
(“r”) period. 


k 

Qk 

Pk 

rk 

0 

1.000000000000000 

- 

- 

1 

4.400000000000000 

35.574 

0.2272727273 

2 

13.43750000000000 

8.783 

0.3274418605 

3 

59.99895833333333 

3.171 

0.2239622216 

4 

318.5420138888889 

1.361 

0.1883549288 

5 

1665.565227141204 

0.575 

0.1912515996 

6 

8563.762388478974 

0.239 

0.1944898926 

7 

43904.44785527987 

0.099 

0.1950545516 

8 

225048.4315368017 

0.041 

0.1950888862 

9 

1153503.664577334 

0.017 

0.1950998843 

10 

5912167.013677382 

0.007 

0.1951067454 

11 

30302311.46598434 

0.003 

0.1951061397 

12 

155313062.3650712 

0.001 

0.1951047195 


TABLE X. Coefficients of the Fade approximants of order 6 of the radial and longitudinal 
perturbations of the variational curve due to the parallactic terms with synodic (“r”) period. 


k 

4 

h" 

n'" 


0 

1.000000000000000 

1.000000000000000 

1.000000000000000 

1.000000000000000 

1 

-2.698270836202914 

-7.098270836202914 

-1.839002476253125 

-7.039002476253125 

2 

-6.296177502041919 

11.49871417725090 

-7.776514253041273 

11.36379862347498 

3 

9.635571285915953 

-5.574715065844693 

4.712486681328506 

-6.420060419171178 

4 

11.54422294537151 

-11.09416027055608 

8.382766821277098 

-10.90642205141127 

5 

-7.427911939834412 

21.91801429175468 

-4.643821407637831 

18.97128636175652 

6 

-10.77769175983123 

-27.61496107523247 

-3.937875173153431 

-20.15607339834716 
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APPENDIX C: CHARACTERISTIC MULTIPLIERS, COMMENSURABILITIES 

AND INSTABILITY 


Let us first recall the basic concept of characteristic multipliers. The small perturbations 
around a periodic orbit in the restricted three-body problem can be described in terms of 
a two-dimensional Poincare map: this is the application connecting two successive intersec¬ 
tions of the trajectory in phase space with a two-plane transversal to the orbit. [One works, 
say, in the three-dimensional reduced phase space corresponding to a hxed value of the Ja¬ 
cobi integral; see for a catalog of such Poincare maps in the case of the Hill problem.] 
For inhnitesimal perturbations, the Poincare map reduces to a linear transformation of the 
plane, leaving hxed the origin which corresponds to the reference periodic orbit. The two 
eigenvalues (Ai, A 2 ) of the inhnitesimal Poincare map (a 2 x 2 matrix) are the characteristic 
multipliers. From the Hamiltonian nature of the dynamics it follows that these multipliers 
are either of the form (e*",e“*") or ,ee~^) with e = ±1 |^. The hrst case means 
generally (apart from the exceptional cases where a = 27r/3, or where a = 27r/4 and some 
inequality is not satished) that the periodic orbit is (quasi-)stable. The second case means 
that the periodic orbit is unstable. A useful quantity for studying the possible loss of sta¬ 
bility is half the sum of the multipliers: a = |(Ai -|- A 2 ), which is either cos a (in which case 
|a| < 1) or £cosh/d (in which case |a| > 1). The loss of stability can occur (apart from the 
above mentioned exceptional cases) only when a crosses the values ±1. The stability of all 
the periodic orbits in Hill’s problem has been studied by Henon We are interested here 
in the families “g'” and of periodic orbits which correspond, respectively, to prograde and 
retrograde lunar-type orbits. Henon found that retrograde orbits are stable for all values of 
the Jacobi integral (i.e. — 1 < m < 0), while prograde orbits are stable only for close enough 
orbits, 0 < m < mcr with|^ 


me,. = 0.1951039966... . (Cl) 

For this value a crosses the value J-l. 

Let us translate this result in terms of the perigee precession of a perturbed Hill or¬ 
bit. Perturbations of Hill’s orbit can be described in terms of the (isoenergetic) normal 
displacement q = {X5Y — Y5X )/This variable satishes “Hill’s equation” 

> (C2) 

where 0(r) = Oq + 2J2T cos2jr is periodic with period vr. The source term a^r) on the 
right-hand side of Eq. (|C^) is zero for free perturbations (i.e. corresponding to adding some 
“eccentricity” to Hill’s “circular-like” variational orbit), and non-zero when one perturbs 


In Henon’s fourth paper |^, he gives the value irmcr = T/2 = 0.61294. The more precise value 
(Cl) was privately communicated to us by Henon, and follows also independently from our results 
in Appendix B (study of the geometric-like growth of various m-series and of the zeroes of Fade 
denominators). 
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Hill’s Hamiltonian -Hhui = (kinetic terms) + Fq + F 2 (e.g. by adding the Fi perturbation 
we are mainly concerned with, or the parallactic terms F^ + ■■■). Perigee precession is 
described by the general homogeneous solution of Eq. ( |C2|) (a = 0). The latter general 
homogeneous solution can be written as a linear combination of complex solutions of the 
form g(r) = Yl,j and of their complex conjugates (we recall that ( = e*’’). On the one 
hand, the quantity c is linked to the usual rate of perigee advance dw/dt (in the non-rotating 
frame) by 0 


dl , , / 1 dw^ 

c = — = (l-|-m) l-— 

dr \ n dt , 


(C3) 


(where I = nt + e — w is the “mean anomaly”). On the other hand, the quantity c is directly 
linked to the characteristic multipliers. Indeed, when r —> r -|- 27r, q{T) gets multiplied by 
g27ric, gg|;g multiplied by Therefore the half-sum of the multipliers is simply 


a = cos(27rc) . 


(04) 


+ 


By using the perturbative series giving the perigee precession, n ^{dzu/dt) = ^ 

-!-■■■ (which has been computed to high accuracy in Ref. [^; see also [0), one can 


check that the crossing of a = -|-1 found by Henon corresponds, when m increases from 
0 to rUcr, to c increasing from 1 to a slightly higher value (~ 1.1) and then decreasing to 
reach the value 1 at m = mcr- From the smoothness of the variation of the characteristic 
multipliers, and therefore of a, with m we deduce that, beyond m = rricr, c{m) goes through 
a quadratic branch point c(m) ~ (1 — m/rricrY^'^ and becomes complex]^. 

Finally, the important information for our purpose is that when m increases up to rricr, the 
quantity c(m) ~ (1 — m/rricrY^'^ is such that both functions cos27rc(m) (= a) and cos7rc(m) 
cross smoothly (without branch points or discontinuities of derivatives) their corresponding 
limiting values cos27r = -|-1 and costt = — 1. 

Let us now consider Hill’s way of solving equation (|C^) . When inserting g(r) = 
into (|C^) one gets an inhnite system of linear equations for the coefficients bj. 
When written in a suitably normalized way, the determinant of this inhnite system (which 
depends on c), say A(c), is a well-dehned quantity (Hill’s determinant). An homogeneous 
solution {a = 0) exists only for the values of c for which A(c) = 0. On the other hand, if we 
consider the case where there is a source term on the right-hand side of Eq. (|C^) of the form 
(j(r) = the corresponding inhomogeneous solution ginhom('r) will have the same 

form as o'(r), but will contain, as for usual hnite determinants, a factor l/A(co-). This fac¬ 
tor will become inhnite when c^ tends to one of the free perigee precession values for which 
A(c) = 0. The analysis of the determinant A(c) = 0 [^, shows that, as a function of c. 


^®The combined facts that rricr is rather small and that c(m) has a quadratic branch point at 
m = trier “explains” the notoriously bad convergence of the perturbation series giving dw/dt. 
Rewriting this series in terms of better behaved quantities, such as cos(27rc) or cos(7rc), improves 
more its convergence than the “Euler transformations”, m m/(l -|- am), which have been 
traditionally used [^. 
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it is a linear function of costtc. We conclude that the forced perturbation ginhom('r) will con¬ 
tain the factor (costtCo- — cos 7 rc(m))“^, where c(m) denotes the free perigee precession value. 
From our analysis above, we know that when m crosses rricr, cos vrc crosses smoothly —1. The 
hnal conclusion is that the source terms for which costtCo- = —1, i.e. Cg- = ±1, ±3, ±5,.. 
generates normal displacements of Hill’s orbit which have pole singularities oc (m — mcr)~^ 
as m crosses the value (|Cl|) found in linear stability analyses. Moreover, it is easily checked 
that the addition of a perturbing potential, say Fp, to Hill’s potential Fq + F 2 , generates a 
source term (T(r) in Hill’s variational equation (|C^) which is a linear combination of Fp and 
3?e [Du{dFp/du)] with real coefficients of the form ko + 2 J2j kj cos2jr. More precisely 

(T(r) = —2(p |(p^<FFp -|- 3?e [Du [dFp/du)]^ , (C5) 

where 

= —DuDu , (C 6 a) 

$(«,«) =mDuDu — 2‘Ste [Dn (^Fniu/Sn)] , (C 6 b) 

with Fniii = Fq + F 2 . 

Therefore, if Fp = will have the form aij) = Ici 

other words, = ±Cp, so that the perturbing potentials Fi, F 3 , F 5 , ... generate source 
terms with Cg- = ±1, Co- = ±3, c^- = ±5, ..., respectively. 


APPENDIX D: LAPLACE ON THE EQUIVALENCE PRINCIPLE 


In the hrst volume of his Traite de Mecanique Celeste s (presented to the French 


Academy of Sciences in 1799), Laplace lists a series of facts suggesting that gravity is pro¬ 
portional to the masses. This list (which is probably inspired by a corresponding list in the 
Principia, although Laplace does not mention Newton here) contains Newton’s argument 
that the motion of satellites would be very sensitive to a violation of the universality of free 
fall, but does not quantify it. As far as we are aware, the only quantitative work of Laplace 
on this idea [^] is contained in the last book of the last volume of the Traite de Mecanique 
Celeste which was presented to the Academie des Sciences on 16 August 1825 . 

The fact that Laplace was then 76 year old (he died a year and a half later on 5 March 
1827) may explain why this work of Laplace contains some strange leaps of reasoning^ It 
seems plausible to us that Laplace, when writing this chapter, was using previous notes of 
his which contained more detailed calculations and more consistent reasonings. Anyway, the 
aim of this Appendix is to show that the hnal limit he quotes on a possible violation of the 
equivalence principle. 


|<^ 12 | < 


3410000 


~ 2.9 X 10 


-7 


(Dl) 


Basically he mentions that his lunar theory and combined (selected) lunar and solar parallax 
data agree to about 1.2 % and then goes on to admit 1/8 = 12.5 % as fractional upper limit on 
the theory/observation agreement. 
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is a very reasonable (and slightly pessimistic) bound, which can be derived in a logically 
clear way using only the information Laplace had in hand. 

Laplace’s new idea (compared to Newton) was to use the “parallactic inequality” in the 
longitude of the Moon as a sensitive test of the equivalence principle[5 “Parallactic inequal¬ 
ity” means the coefficient of the synodic term — sin r in the expression of the lunar longitude 
n as a function of time[^. We have computed this coefficient, say A, in Appendices A and 
B for the Main Problem (i.e. neglecting eccentricities and inclinations). Its theoretical 
expression reads 

= ApAR -l- Arp , (D2) 

where the normal “parallactic” contribution reads 

ApAR = -^{^2 — Wi) S'pAR(m) , (D3) 

8 a' 

while the “equivalence principle” one reads 


Arp — 


3m 5 i2 — *S'Rp(m) . 
a 


(D4) 


Here, S'pAR(m) = 1 -|- ^m S'Rp(m) = 1 -|- ym are slowly converging series in 

powers of m = n'/{n — n'). In the third volume of his Traite de Mecanique Celeste, Laplace 
computes Arar with particular care, pushing the calculation to hfth order in m inclusively^^. 
He was therefore entitled to considering that the theoretical error on Arar was negligible 
compared to the observational uncertainties in A°*”^. Note that Arar is proportional to 
the inverse distance to the Sun, i.e. to the “solar parallax” ns = Re/o.', where Re is the 
equatorial radius of the Earth (hence the name “parallactic inequality”). The result of 
Laplace can be expressed as 


Arar — 14.3 ns , 


(D5) 


where both Arar and ns are expressed in seconds of arc[^ 

On the observational side, Laplace had in hand both lunar data and data on the solar 
parallax. The two phenomenological “lunar tables” he was using (one by Mason and one by 


Note that as early as 1753, T. Mayer had used the theory of the parallactic inequality to infer 
the value of the solar parallax . 

Beware that Laplace was actually working with the inverse function: t = t(v). 

We have checked the first orders of his result and found them to agree with ours. Note that 
Laplace includes the effect of eccentricities and inclinations that we neglect. 

Beware that Laplace, in his volume III, uses (Revolution oblige) “decimal seconds”, i.e. 10“® 
of a right angle. 
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Burg) gave for the “observed” value of the synodic inequality in longitude, ^Mason = 116.’ 
and = 122.” 38 


33 . As for the solar parallax, many scientihc expeditions had taken 


advantage of the passages of Venus in front of the Sun in 1761 and, especially, 1769 to measure 

Q for a more exhaustive 


TTs (see 1^ for a detailed text by a contemporary of Laplace, or 
historical treatment). The published results ranged between = 8.”43 (Planman) and 
7 r°’’® = 8.”80 (Pingre) 0, The comparison between the theoretical results (p^ , (pSf) 
and the observational results on A and gives a value for a possible equivalence-principle- 
violation contribution 


Aep = - 14.37r°^" . (D6) 

Worst-case limits on App are obtained by taking the extreme values on the right-hand 
side of (P^) (e.g. = Ag’jfj.g — 14.37rg'®'“™'^“)Q. This yields 

-9.”2 < Aep < l.”8 . (D7) 

On the other hand, the theoretical result (|D4|) [using a'/a ~ 391 and S'Ep(m) ~ 1.72 (see 
Appendix B)] reads Aep — 163 5i2 — 3.36 (5i2 x 10^ in seconds of arc, so that we get the 
following worst-case bounds on 612 

- 2.7 X 10“^ < ^12 < 0.54 X 10“^ . (D8) 

From this point of view, the hnal bound quoted by Laplace, Eq. (0), seems very reasonable 
and consistent with the observational uncertainties in his time. Note, however, that Laplace 
never quotes a precise theoretical formula for Aep. He only says (and uses) the fact that 
the synodic amplitude A is proportional to its source term in the perturbing function. This 
neglects the leading “p dependence” of oc 2-|-p in Eq. (A19b|) which says that Aep/Apar x 
3/5 X S'Ep(m)/S'pAR(r«). In Laplace’s published analysis the lacking (unfavourable) factor 
3/5 is effectively compensated by his overpessimistic estimate of the fractional uncertainty 
on A: 1/8 = 12.5 Vcf]. 

Let us end this section by raising a historical question. Though Laplace was fully aware of 
the scientihc interest of the bound ( |D1|) , and of the fact that it was (at the time) more precise 
than the bounds obtained from ground tests of the universality of free fall, his successors in 
celestial mechanics seem (as far as we know) to have lost interest in the issue. However, near 
the end of the nineteenth century, especially after the theoretical work of Hansen, Delaunay, 
Hill and Brown (who improved the computation of Apar), and after many improvements 
on the observational side, it should have been possible to obtain more stringent limits on 
di 2 - For instance, Delaunay computed A^p to which corresponds to a truncature error 


Such a worst-case approach seems appropriate to a pre-least-squares-law period. Before Gauss’ 
theory of measurement errors, scientists quoted only “central values” for measured quantities. 

As for the subleading dependence on p, i.e. the ratio Sep (itt.)/Spar( m) ~ 1.72/1.60 ~ 1.08 
Laplace’s experience with similar factors in many terms of lunar theory might have suggested him 
that he did not need to worry about it. 
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of 0.”3366/125.”4201 ~ 2.7 x 10“^ (for the principal part of A; see also Appendix B). 
This is negligible (when added in quadrature) compared to the observational error on vr^ 
at the time. For instance, the laboratory measurements of the velocity of light by Foucault 
and Cornu gave values ranging between = 8.”834 and tts = 8.”881, the passage of 
Venus across the solar disk in 1874 gave a range 8.” 76 — 8.” 88 and the recommended 
value starting in 1896 was 8.” 80. This suggests that a reasonable upper bound on the 
uncertainty of tts at that time was 0.”08 (i.e. 0.9 %). [By comparison, the modern value is 
TVs = 8.”794148 [^.] On the other hand, the observational error on at the time was 
< I(125.”46 — 124.” 70) = 0.”38, i.e. < 0.3 % (see [^, p. 533). This has a negligible effect 
on the derivation of a bound on Si]^. Using Eq. ( |D6| ) with A^p — 3.36 5i2 x 10^ seconds of 
arc, this leads to the bound |Ji 2 | < 14.3(57rg'^'^/(3.36 x 10’^) ~ 3.4 x 10“®, which is slightly 
better than the value obtained by Eotvos in 1890 (I^abI < 5 x 10“® []^). We do not wish 
to take too seriously such a posteriori derivations of limits on 612 , but we consider this as 
an interesting example of historical eclipse of a deep concept, which has been rejuvenated, 
within a new theoretical and observational context, only in the last decadesf^. 


Let us note for completeness that the determinations of the mass of the Moon at the time were 
accurate enough to estimate with negligible error the mass ratio factor X 2 — Xi entering Apar. 
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For completeness, let us note that Poincare kept alive 


result in one of his popular books |64|. 


this concept by mentioning Laplace’s 


47 







REFERENCES 


* On leave from the Institute of Astronomy, Charles University, Svedska 8, 15000 Prague 
5, Czech Republic (e-mail: davok@aci.cvut.cz). 

[1] G. Galilei, Discorsi e dimostrazione matematiche intorno d due nuove scienze, (Elzevir, 
Leiden, 1638); French translation by M. Clavelin, Discours concernant deux sciences 
nouvelles (A. Cohn, Paris, 1970); see p. 61: “... it came to my mind that, if the resistance 
from the medium were totally suppressed, all bodies would go down with the same 
speed.”; see also p. 70 where Galileo claims to have verihed the universality of free fall, 
to better than 0.1 %, by pendulum experiments. The pages quoted here correspond, 
respectively, to p.ll6 and p.l29 of the Edizione nazionale, and to p.72 and p.85 of the 
English translation by H. Crew and A. de Salvio, Dialogues concerning two new sciences, 
(Dover, New York, 1954). 

[2] Though this issue has a long prehistory (starting, say, with J. Philoponos, 5th Century 
A.D.), even G.B. Benedetti who came closest to stating clearly the universality of free 
fall (in his Resolutio of 1553, and his Demonstratio of 1554) did not quite succeed in 
doing so. See Mechanics in Sixteenth-Century Italy; Selections from Tartaglia, Benedetti, 
Guido Ubaldo and Galileo, translated and annotated by S. Drake and I.E. Drabkin (The 
University of Wisconsin Press, Madison, 1969). 

[3] I. Newton, Philosophiae Naturalis Principia Mathematica (London, 1687); English 
translation by F. Cajori, Newton’s Principia (University of California Press, Berkeley, 
1934). 

[4] See Book III, Proposition VI, Theorem VI of the Principia Mathematica, pp. 411- 
413 of Cajori’s edition; see also the discussion in T. Damour, chapter 6 of 300 years 
of gravitation, edited by S.W. Hawking and W. Israel (Cambridge University Press, 
Cambridge, 1987), pp. 128-198, notably the Section 6.6, pp. 141-144. 

[5] P.S. Laplace, Traite de Mecanique Celeste, Volume V (Bachelier, Paris, 1825). Laplace’s 
work is contained in Book XVI, Chapter IV. Reprinted in Oeuvres de Laplace, Tome 
Cinquieme (Gauthier-Villars, Paris, 1882), pp. 445-452. 

[6] F.W. Bessel, Ann. der Physik und Chemie (Poggendorff) 25, 401 (1832). 

[7] R. V. Eotvos, Math. u. Naturw. Ber. aus Ungarn 8, 65 (1890). 

[8] R. V. Eotvos, D. Pekar and E. Fekete, Ann. der Physik 68, 11 (1922). 

[9] A. Einstein, Jahrb. Radioact. Elekt. 4, 411 (1907). 

[10] K. Nordtvedt, Phys. Rev. 170, 1186 (1968). 

[11] K. Nordtvedt, Phys. Rev. 169, 1014 (1968). 

[12] K. Nordtvedt, Phys. Rev. 169, 1017 (1968). 

[13] R.H. Dicke, Rev. Mod. Phys. 29, 355 (1957); see notably pp. 356, 357 in which by “weak 
interaction energy” Dicke means both Fermi interactions and gravitation. 

[14] R.H. Dicke, in Evidence for Gravitational Theories, Proceedings of the International 
School of Physics “Enrico Fermi”, Course XX, edited by C. Mqller (Academic Press, 
New York, 1962), pp. 1-49, see notably pp. 16-19 and 26-27. 

[15] R.H. Dicke, in Gravitation and Relativity, edited by H.Y. Chiu and W.F. Hoffman 
(Benjamin, New York, 1964), pp. 1-16, see notably pp. 13-16. 

[16] K. Nordtvedt, Phys. Rev. D 7, 2347 (1973). 


48 



[17] D. Brouwer and G.M. Clemence, Methods of Celestial Mechanics, (Academic Press, 
Orlando, 1961). 

[18] V.A. Brumberg and T.V. Ivanova, Trans. Inst. Theor. Astron. (Leningrad) 19, 3 (1985). 

[19] C.M. Will, Theory and Experiment in Gravitational Physics, (Cambridge University 
Press, Cambridge, 1981); second edition 1993. 

[20] J.C. Williams et al., Phys. Rev. Lett. 36, 551 (1976). 

[21] LI. Shapiro et ah, Phys. Rev. Lett. 36, 555 (1976). 

[22] J.O. Dickey, X.X. Newhall and J.C. Williams, Adv. Space Res. 9, No. 9, 75 (1989). 

[23] J. Muller, M. Schneider, M. Soffel and H. Ruder, Astroph. J. 382, LlOl (1991). 

[24] C. de Pontecoulant, Theorie analytique du systeme du monde. Volume IV, (Paris, 1846) 

[25] K. Nordtvedt, Icarus 114, 51 (1995). 

[26] W. de Sitter, Mon. Not. R. Astr. Soc. 76, 699; 77, 155 (1916). 

[27] V.A. Brumberg, Bull. Inst. Theor. Astron. (Leningrad, USSR) 6, 733 (1958). 

[28] M.H. Soffel, Relativity in Astrometry, Celestial Mechanics and Geodesy, (Springer, 
Berlin, 1989). 

[29] V.A. Brumberg, Essential Relativistic Celestial Mechanics, (Adam Hilger, Bristol, 
1991). 

[30] J.F. Lestrade and M. Chapront-Touze, Astron. Astroph. 116, 75 (1982). 

[31] M. Soffel, H. Ruder and M. Schneider, Astron. Astroph. 157, 357 (1986). 

[32] B. Shahid-Saless, Phys. Rev. D 46, 5404 (1992). 

[33] J.C. Williams, X.X. Newhall and J.O. Dickey, Phys. Rev. D, in press. 

[34] T. Damour and G. Esposito-Farese, Phys. Rev. D 50, 2381 (1994). 

[35] T. Damour, in Gravitation and Quantizations, Proceedings of the LVIIth Les Houches 
Summer School, edited by B. Julia and J. Zinn-Justin (North-Holland, Amsterdam, 
1995), pp. 1-62 . 

[36] T. Damour and A.M. Polyakov, Nucl. Phys. B 423, 532 (1994); and Gen. Rel. Grav. 
26, 1171 (1994). 

[37] T. Damour and G. Esposito-Farese, Class. Quantum Grav. 9, 2093 (1992). 

[38] T. Damour, M. Soffel and C. Xu, Phys. Rev. D 43, 3273 (1991); 45, 1017 (1992); 47, 
3124 (1993); 49, 618 (1994); see also T. Damour and D. Vokrouhlicky, submitted to 
Phys. Rev. D. 

[39] G.W. Hill, Amer. J. Math. 1, 5 (1878); 1, 129 (1878); 1, 245 (1878); Acta Math. 8, 1 
(1886). 

[40] E.W. Brown, An Introductory Treatise on the Lunar Theory, (Cambridge University 
Press, Cambridge, 1896); reprinted by (Dover, New York, 1960). 

[41] A. Wintner, Math. Z. 24, 259 (1926). 

[42] Y. Hagihara, Celestial Mechanics, Volume II, Part 2 (The MIT Press, Cambridge, 1972). 

[43] C.L. Siegel and J.K. Moser, Lectures on Celestial Mechanics, (Springer-Verlag, Berlin, 
1971). 

[44] A.M. Liapunov, Trudi Otd. Fys. Nauk Ljobit. Estestvoznan. 8, 1 (1896); reprinted in 
Liapunov Sobranie Sochinenii I, 418 (Akad. Nauk SSSR, 1954). 

[45] H. Poincare, Legons de Mecanique Celeste, Volume 2, Part 2, (Gauthier-Villars, Paris, 
1909). 

[46] M. Moons, MINIMS: User’s guide, (Rapport interne No. 91/15 du Departement de 
Mathematique FUNDP, Namur, 1991). 


49 



[47] M.C. Gutzwiller and D.S. Schmidt, Astronomical Papers prepared for the use of the 
American Ephemeris and Nautical Almanach, Volume XXIII, Part I (1986). 

[48] Connaissance des Temps pour I’An 1995, Bureau des Longitudes, (Service Hydro- 
graphique et Oceanographique de la Marine, Brest , 1994). 

[49] M. Henon, Ann. Astr. 28, 499 (1965); 28, 992 (1965); Bull. Astr. Paris 1, fasc. 1, 57 
(1966); 1, fasc. 2, 49 (1966); Astron. Astroph. 1, 223 (1969); Astron. Astroph. 9, 24 
(1970). Hill’s problem is treated in the last two references. 

[50] J.O. Dickey et ah. Science 265, 482 (1994). 

[51] J. Laskar and O. Neron de Surgy, private communication (1995). 

[52] Y. Su et ah, Phys. Rev D 50, 3614 (1994). 

[53] T. Damour and K. Nordtvedt, Phys. Rev. Lett. 70, 2217 (1993); Phys. Rev. D 48, 3436 
(1993). 

[54] B. Chauvineau and F. Mignard, Atlas of the Circular Planar Hill’s problem, (Internal 
Publication of the Observatoire de la Cote d’Azur, Grasse, 1991). 

[55] A. Deprit, J. Henrard and A. Rom, Astron. J. 76, 273 (1971). 

[56] J. Henrard, Celest. Mechanics 17, 195 (1978); Ciel et Terre, 1-27 (1973). 

[57] V.I. Arnold, Mathematical Methods of Classical Mechanics, (Springer-Verlag, New York, 
1989). 

[58] P.S. Laplace, Traite de Mecanique Celeste, Volume I (Paris, 1799), hrst part. Book 
H. Reprinted in Oeuvres de Laplace, Tome Premier (Imprimerie Royale, Paris, 1843), 
p. 139. 

[59] C.C. Gillipsie, editor. Dictionary of Scientific Bibliography, Volume XV, Supplement I 
(Scribner, New York, 1970), pp. 387-403. 

[60] P.S. Laplace, Traite de Mecanique Celeste, Volume HI (Paris, 1802), Book VH. 
Reprinted in Oeuvres de Laplace, Tome Troisieme (Imprimerie Royale, Paris, 1844), 
pp. 196-327. 

[61] J.C. Houzeau, Vace-Mecum de TAstronome, Annales de 1’Observatoire Royal de Brux¬ 
elles, (F. Hayez, Bruxelles, 1882), pp. 403-411. 

[62] J. Lalande, Astronomie, Tome Seeond, (Didot I’Aine, Paris, 1792). 

[63] H. Woolf, The transits of Venus, (Princeton University Press, Princeton, 1959). 

[64] H. Poincare, Science et Methode (Flammarion, Paris, 1908); see book HI, chapter HI, 
section HI. 


50 




d 


0.8 - 0.6 - 0.4 - 0.2 0 0.2 






